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'd\'. 

y^ ' Some formulae for the perturbations of the matter fields are summarized within the 

>— ^ , framework of the second-order gauge-invariant cosmological perturbation theory in a four 

^^ ■ dimensional homogeneous isotropic universe, which is developed in the papers [K. Nakamura, 

r^ ' Prog. Theor. Phys. 117 (2007), 17.]. We derive the formulae for the perturbations of the 

C^ ' energy momentum tensors and equations of motion for a perfect fluid, an imperfect fluid, 

^~~i , and a signle scalar field, and show that all equations are derived in terms of gauge-invariant 

pf-\ ■ variables without any gauge fixing. 

'TT'. §1. Introduction 

5^ ' The general relativistic second-order cosmological perturbation theory is one of 

f^. topical subjects in the recent cosmology. By the recent observation of the Cosmic 

\ Microwave Background (CMB) by Wilkinson Microwave Anisotropy Probe, ^^ the 

■^ ■ first order approximation of our universe from a homogeneous isotropic one was re- 

JC . vealed. The observational results suggest that the fluctuations of our universe are 

Tjj- \ adiabatic and Gaussian at least in the first order approximation. As a next step, 

OO ' the clarifications of the accuracy of these observational results are actively discussed 

^ : . both in the observational^ and the theoretical side^^'^ through the non-Gaussianity, 

^T I the non-adiabaticity, and so on. With the increase of precision of the CMB data, 

QQ ■ the study of relativistic cosmological perturbations beyond linear order is a topical 

f^ . subject especially to study the generation of the primordial non-Gaussianity in in- 

^ \ flationary scenarios^' and the non-Gaussian component in CMB anisotropy. ^^ The 

second-order cosmological perturbation theory is one of such perturbation theories 
beyond linear order. 
C^ \ According to this physical motivation, we proposed a clear gauge-invariant for- 

mulation of the second-order general relativistic cosmological perturbation theory.^' 
In this paper, we refer these works as KN2007. This gauge-invariant formulation of 
the second-order cosmological perturbations is a natural extension of the first-order 
gauge-invariant cosmological perturbation theory.^-*"^) The formulation in KN2007 
is one of the applications of the gauge-invariant formulation of the second-order per- 
turbation theory on the generic background spacetime developed in the papers by 
the present author.^)' ^^^ These papers are referred in this paper as KN2003^^ and 
KN2005.^^) This general formulation is a by-product of the investigations of the 
oscillatory behaviors of self-gravitating Nambu-Goto membranes. ^^^ 

In KN2007, we defined the complete set of the gauge-invariant variables of the 
second-order cosmological perturbations in the Friedmann-Robertson- Walker uni- 
verse based on the formulation developed in the papers KN2003 and KN2005. We 
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considered the two cases of the Friedmann-Robertson- Walker universe: one is the 
universe filled with the single perfect fluid and another is the universe filled with the 
single scalar field. We also derived the second-order Einstein equations of cosmo- 
logical perturbations in terms of these gauge-invariant variables without any gauge 
fixing in these two cases. We have also found that the procedure to find gauge invari- 
ant variables proposed in KN2003 plays a crucial role in the derivations. Rather, we 
can use the formulae proposed in KN2003 to check whether the resulting formulae 
are correct or not. 

This paper is the second part of KN2007. In this paper, we summarize the 
formulae for the components of the first- and the second-order perturbations of the 
energy momentum tensors and the equations of motion which are derived from the 
divergence of the energy momentum tensors. As the matter contents, we consider 
the three matter fields: a single perfect fluid; a single imperfect fluid which includes 
additional terms of the energy flux and anisotropic stress to the perfect fluid; and a 
single scalar fleld. In the early universe, photon and neutrinos should be described 
by the Boltzmann distribution functions. '^^'^'^^ Photon's interaction with baryon 
and the free streaming of neutrinos lead anisotropic stress and these effects will be 
important in the recent cosmology. Although the energy flux and the anisotropic 
stress in the imperfect fluid are determined through these micro-physical processes, 
we just phenomenologically treat these terms in this paper. 

Although the perturbative expressions of the energy momentum tensor and equa- 
tions of motion were also derived in some literatures, ^^^ in this paper, we show alter- 
native derivations of these perturbations. In our derivations, we respect the gauge 
invariance of the perturbative variables. We again show that the formulae of gauge 
invariant-variables proposed in KN2003 [Eqs. (j2-20p and (|2-2ip in this paper] also 
play crucial roles in the derivations of perturbative expressions of the equations for 
matter flelds. The first- and the second-order perturbations of the equations of mo- 
tion are decomposed into gauge-invariant and gauge- variant parts as Eqs. (|2-20p and 
(|2-2ip . respectively. In these derivations, we do not fix any gauge degree of freedom. 
In spite of this no gauge-fixing, we show all perturbations of the equations of mo- 
tion are given in gauge-invariant forms through the lower order perturbations of the 
equations of motion for matter fields. In this sense, we may say that the general 
framework of the second-order gauge invariant perturbations proposed in KN2003 
and KN2005 does work not only in the perturbations of the Einstein equations but 
also in the equations of motion for the matter fields. The main purpose of this paper 
is to show this. 

Further, in this paper, we do not ignore the first-order vector- and tensor-modes 
which are ignored in KN2007. Moreover, in the derivation of the perturbations 
of the energy momentum tensors and the equations of motion, we do not use any 
information of the Einstein equations. Therefore, the formulae derived in this paper 
are valid even if we consider any other theory of gravity than the Einstein theory. 

The organization of this paper is as follows. In ^ we briefiy review the defi- 
nitions of the gauge-invariant variables for the first- and second-order perturbations 
which were defined by KN2007.^^ In ^ we derive the first- and the second-order 
perturbations of the energy momentum tensors and equations of motion for a perfect 
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fluid, an imperfect fluid, and a scalar field. In the derivation in this section, we do 
not specify the background spacetime. Therefore, the formulae summarized in this 
section are valid in perturbation theories on any background spacetime. In ^ we 
derive the explicit expression of the components of the energy momentum tensors 
and the equations of motion for matter fields. The final section, ^ is devoted to 
the summary and the discussions. Further, in Appendix [Aj we explicitly give the 
components of the perturbations of the acceleration, expansion, shear, and rotation 
associated with the fluid four-velocity, which are necessary to derive the results in 
§g3]and[l 

We employ the notation of our previous papers KN2003,^'* KN2005,^'^^ and 
KN2007^' and use the abstract index notation. ^^^ We also employ the natural unit 
in which the light velocity is denoted by c = 1. 

§2. Gauge-invariant variables in the second-order perturbations 

In this section, we briefly review the definitions of the gauge-invariant variables. 
First, in N2.ll we review the gauge degree of freedom and the first- and the second- 
order gauge transformation rules. Then, we briefly explain the gauge-invariant vari- 
ables for the metric (in N2.2p and matter perturbations (in N2.3p . 



2.1. Gauge degree of freedom 

In any perturbation theory, we always treat two spacetime manifolds. One is the 
physical spacetime ^A = M.\ and the other is the background spacetime Mq. The 
physical spacetime M.\ \s our nature itself which we want to describe through the 
perturbations. On the other hand, the background spacetime Mq is just a reference 
spacetime for the calculations of perturbations. Although this background spacetime 
has nothing to do with our nature, to calculate perturbations, it is necessary to 
introduce this reference spacetime Mq by hand. Since these two spacetime manifolds 
are distinct from each other, we have to introduce a point-identification map X\ : 
Aio -^ M.\. This point-identification map Xx called a gauge choice in perturbation 
theories. Through the pull-back X^ of the gauge choice Xx, any physical variable Q\ 
on the physical manifold M.\ \s pulled back to X^Qx on the background spacetime 
A^Q. The pull-back X^Qx is a representation on the background spacetime M.q of 
the physical variable Qx on the physical spacetime Mx- Although we do not know 
about the physical spacetime M.x ^^ the starting point of the perturbation theory, 
we can treat the physical variable Qx on the physical manifold Mx as the variable 
'^aOa on the background spacetime A^o through this pull-back X^. 

In the case of the perturbations in the theory with general covariance, the above 
strategy of the perturbation theory includes an important trouble. This is the fact 
that the gauge choice Xx is not unique by virtue of general covariance. Rather, 
there is degree of freedom in the gauge choice Xx, i.e., we may apply the different 
point-identification map yx from Xx as a gauge choice. In this case, the represen- 
tation yxQx on the background spacetime Mq of the physical variable Qx on the 
physical spacetime A^a is different from the representation X^Qx. This difference is 
unphysical because it has nothing to do with the nature of the physical spacetime 
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Ai\. We can also consider the transformation rule from a gauge choice A\ to another 
one y\, which is called gauge transformation. The gauge transformation X\ -^ yx 
is induced by the diffeomorphism ^x '■= (-^a)^^ ° 3^A- Actually, the diffeomorphism 
'Px changes the point-identification maps from Xx to 3^;^. Further, the pull-back <^^ 
of the diffeomorphism ^x changes the representation X^Qx of the physical variable 
Qx to another representation y^Qx'- 

y*xQx = y*x i^x ° '^'a^')* Oa = yiiXx'r^xQx = n^iQx. (2-1) 

Since ^x is the diffeomorphism on the background manifold Mq^ the Taylor 
expansion of the pull-back ^^ is given by 

^a-^a^Oa = -^a*Oa + A^A'^Qa + \\^ (% + £l) X*xQx + 0(A3), (2-2) 

where 5,1 and ^2 ^^^ generators of the diffeomorphism c^^.^^) On the other hand, we 
consider the perturbative expansion 

<^a*Qa = Qo + aJ^Q + ]^\^fQ + 0(a3) (2-3) 

of the representation X^Qx under the gauge choice Xx^ where Qo is the background 
value of the variable Qx- The first- and the second-order perturbations -^Q and 
^"^'xQ are defined by this equation (j2-3p . From Eqs. (|2-2p and ()2-3p . we can easily 
derive the gauge transformation rule of each order: 

^i^Q-^5Q = %)Qo, (2-4) 

(Jq - (Jq = 2£^J'^Q + {£5,,, + 4,J go. (2-5) 

Further, we introduce the concept of "i/ie order by order gauge invariance" . We call 

(p) 
the pth-order perturbation ^Q is gauge invariant iff 

%^Q = ^$Q (2-6) 

for any gauge choice Xx and yx- We have been considering the concept of "the 
gauge invariant up to order n" in the series of the papers KN2003,^^ KN2005,^'^^ and 
KN2007,^) following the idea by Bruni and Sonego.^^-* However, we should regard the 
gauge invariance in this series of the papers is this "order by order gauge invariance" 
rather than "the gauge invariance up to order n" . This notion of the order by order 
gauge invariance is weaker than the notion of the gauge invariance up to n, since 
we do not say anything about the gauge invariance of the other orders in the above 
order by order gauge invariance. 

Employing the idea of this order by order gauge invariance, we proposed a 
procedure to construct gauge-invariant variables of higher-order perturbations in 
KN2003.^'' Inspecting the gauge transformation rules (j2-4p and (j2-5p . we can de- 
fine the gauge-invariant variables for a metric perturbation and for arbitrary matter 
fields. 
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2.2. Gauge-invariant variables for metric perturbations 

Following the expansion form (I2-3J) . we also expand the pulled-back X^dab of the 
metric gab on the physical spacetime Aiy. 

X^gab = gab + ^Xhah + -^^P^ab + 0{X^), (2-7) 

where gab is the metric on the background spacetime Mq. Although the expansion 
(|2-7p of the metric depends entirely on the gauge choice X\, henceforth, we do 
not explicitly express the index of the gauge choice X\ in expressions if there is 
no possibility of confusion. As shown in KN2007,^^ at least in the cosmological 
perturbation case, the first-order metric perturbation hab is decomposed as 

hab =■ 'Hab + £xgab, (2-8) 

where Tiab and X" are the gauge-invariant and gauge-variant parts of the linear- 
order metric perturbations,^^ i.e., under the gauge transformation (|2-4p . these are 
transformed as 

yHab - xHab = 0, y^"" - xX"" = ^'^^y (2-9) 

Further, the second-order metric perturbation lab is also decomposed as 

lab ='■ Cab + '^£xhab + {^Y - £ x) gab, (2-10) 

where Cab and Y'^ are the gauge-invariant and gauge-variant parts of the second- 
order metric perturbations, i.e., these are transformed as 

yCab-xCab = 0, yY''-xY'' = ^f^^ + [^ii),X]'' (2-11) 

under the gauge transformation ()2-5p . 

In KN2007,^^ the details of the derivation of this gauge-invariant part of the 
second-order metric perturbation are explained in the context of cosmological pertur- 
bations. In the case of the cosmological perturbations, we consider the homogeneous 
isotropic background spacetime whose metric is given by 

gab = a^ {-idr])aidr])b + J^j{dx')a{dx^)b] , (2-12) 

where jab •= lij{dx^)a{dx^)b is the metric on the maximally symmetric three space. 
As shown in KN2007, the decomposition (j2-8p is accomplished if we assume the 
existence of the Green functions Z\"^ := (i:>*A)"\ (^ + 2A')"\ and (Z\ + 3ir)"\ 
where Di is the covariant derivative associated with the metric 7jj on the maximally 
symmetric three space and K is the curvature constant of this maximally symmetric 
three space. We also showed in KN2007 that we may choose the components of the 
gauge-invariant part Tiab of the first-order metric perturbation as 

2 r . (y . , . . , . , . (1) 



Tiab = a^\-2 <p {dr])aidr])b + 2 v^ {dr,)^a{dx')b) 

( (1) (1) \ -1 

+ ( -2 1^ 7i,+ X^A (dx')a{dx3)by (2-13) 
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(1) (1) 

where i' i and Xij satisfy the properties 

D^^l\:=fW,^l\=0, x/:= 7^4V,= 0, D'X^,= 0. (2-14) 

Further, we may also choose the components of the gauge-invariant part Cab of the 
second-order metric perturbation as 

2 f . (2) _ _ _ _ (2) 



Cab = a <^ -2 <? {d'n)a{dv)b + 2 Vi {d'n)(^a{dx')b) 

( (2) (2) \ .1 

+ ( -2 q, 7,j+ XiA {dx')a{dx^)by (2-15) 

(2) (2) 

where v i and Xij satisfy the properties 

• (2) ^^''. ■ (2) 

D'yi=Q, x/=0, D'xi3=^- (2-16) 

Here, we also note the fact that the definitions (j2-8p and (12-lOp of the gauge- 
invariant variables are not unique. This comes from the fact that we can always 
construct new gauge-invariant quantities by the combination of the gauge-invariant 

(1) (1) 

variables. For example, using the gauge-invariant variables ^ and v^ of the first- 
order metric perturbation, we can define a vector field Za by 

Za:= a <^ {dri)a + a vi {dx')a. (2-17) 

Although there is no specific physical meaning in this vector field Za, at least, we 
can say that the vector field Za defined by (|2T7p is gauge-invariant. Through this 
gauge-invariant vector field Za, we can rewrite the decomposition formula (j2-8p for 
the linear-order metric perturbation as 

hab = Ti-ab — £zgab + £zgab + ^XQab, 

= '■ K^ab + £x+Zgab, (2-18) 

where we have defined new gauge-invariant variable )Cab by 

A^afe := "Hab - £zgab- (2T9) 

Clearly, /C^fe is gauge-invariant and the vector field X"" + Z"" satisfies the gauge trans- 
formation rule ()2-9p for the gauge- variant part of the first-order metric perturbations. 
Although the definition of the gauge-invariant variables is not unique, we can specify 
the components of the gauge-variant part Xa without ambiguities if we specify the 
components of the gauge-invariant part TCab as shown in KN2007. In this paper, we 
specify the components of the tensor Hab as Eq. (j2-13p . which is the gauge-invariant 
part of the linear-order metric perturbation associated with the longitudinal gauge. 
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2.3. Gauge-invariant variables for matter fields 

As shown in KN2003, using the above first- and second-order gauge-variant 
parts, X"" and Y", of the metric perturbations, we can define the gauge-invariant 
variables for an arbitrary field Q other than the metric. These definitions imply that 
the first- and the second-order perturbations ^^'Q and ^"^'Q are always decomposed 
into gauge-invariant part and gauge-variant part as 

WQ=:a)Q+£^Qo, (2-20) 

(2)q =: (2)g + 2£x'^^^Q + {£y - £\] Qo, (2-21) 

respectively, where ^^'Q and ^"^'Q are gauge-invariant parts of the first- and the second- 
order perturbations of ^'^'Q and ^"^'Q, respectively. 

Through the formulae (j2-20p and (j2-2ip . we can define the gauge-invariant vari- 
ables for the matter field. In this paper, we consider the cases of a perfect fluid; 
an imperfect fluid; and a scalar field. All these matter fields consist of fundamental 
quantities. For example, we regard the energy density, the pressure, and the four- 
velocity as fundamental variables for a perfect fluid. Here, we show the definitions 
of the gauge-invariant variables for these fundamental quantities. These definitions 
are just following to the formulae (j2-20p and (|2-2ip . However, based on these defini- 
tions of gauge-invariant variables, we will show that all perturbative quantities are 
decomposed into gauge-invariant and gauge- variant parts as Eqs. (|2-20p and (|2-2ip . 

2.3.1. Perfect fluid 

Here, we consider the definitions of the gauge-invariant variables for the funda- 
mental variables of a perfect fluid. The total energy momentum tenor of the fluid is 
given by 

ipyTj = {-e + p)uau'' + p6j, (2-22) 

where the fundamental variables for a perfect fluid are the energy density e, the 
pressure p, and the four- velocity u". We expand these fundamental variables as 

(1) 1 o (2) o 

e~e + Xe +-X^e +0{X^), (2-23) 

p~p + X%\lx^%\o{X^), (2-24) 

Ua -Ua + X (Ua) +-X^ (uj +0{X^), (2-25) 

where e, p, and Ua characterize the pull-back of the fluid on the physical spacetime to 
the background spacetime through an appropriate gauge choice Xx, while e, p, and Ua 
are their background values on the background spacetime. Following to Eqs. (j2-20p 
and (j2-2ip . we define the gauge-invariant variable for the perturbation of the fluid 
components e, p, and Ua'. 

(1) (1) (1) (1) (1) (1) 

S:=e -£xe, V.= P -£xP, Ua:=iua) -£xUa, (2-26) 

fjf -2£x ^6^ - {£y - £'x} e, V:=P -2£x P - {£y - £'x} P, (2-27) 
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(2) (2) (1) ^ 

Ua: = iUa)-2£xUa-{£Y-£x}ua- (2-28) 

2.3.2. Imperfect fluid 

Here, we consider the generic case of an imperfect fluid. The energy-momentum 
tensor is decomposed into fluid quantities based on the orthogonality to the four- 
vector field u"" as 

T; ^ = eUaU^ + P (<5^^ + UaU^^ + qaU^ + UaQ^ + ttJ (2-29) 

= (P>f, ^ + g^'qaUc + g'^uaqc + g'^Ttac, (2-30) 

where 

U^Qa = 0, (2-31) 

7f[,6]=0, 72"7f,fe = 0, 7f/ = ^''Vfafe = 0. (2-32) 

The energy density e, the isotropic pressure p, and the four- velocity Ua of the imper- 
fect fluid in Eq. (j2-29p are expanded as Eqs. ()2-23p ~ (j2-25p and the gauge-invariant 
variables for their perturbations are defined by Eqs. (|2-25p - (|2-28p as in ^2.3.1[ In 
addition to these fluid components, in the imperfect fluid case, we add the energy 
flux Qa and the anisotropic stress Ttab associated with the vector field Ua as fiuid com- 
ponents. Since the first two terms in Eq. (j2-29p coincide with the energy momentum 
tensor for a perfect fluid, we call these terms as "i/ie perfect parf and denote them by 
(pyfji g^g jjj gq_ (j2-3Up . On the other hand, we call the remaining terms in Eq. (|2-30p 

(*>f/ := g,u^ + n,g-^ + 7f,^ (2-33) 

as "i/ie imperfect parf of the energy momentum tensor for an imperfect fiuid. 

Now, we consider the perturbative expansion of the energy fiux g^ and the 
anisotropic stress itab- Although these quantities should be given through the micro- 
physical process, in this paper, we regard these variables as fundamental quantities 
for an imperfect fiuid and expand these variables as 

(1) 1 , (2) 

qa =:qa + X (qa) +-X^ (qa) +0(A=^), (2-34) 

(1) 1 (2) 

TTab =■■ TTab + A (TTab) +-X^ {iTab) +0{X^). (2-35) 

Further, we introduce gauge-invariant variables for the perturbations of the energy 
flux Qa and the anisotropic stress iTab- Following to the decompositions (|2-20p and 
(j2-2ip for an arbitrary matter fleld, the first- and the second-order perturbations of 
the energy fiux and the anisotropic stress are decomposed into gauge-invariant and 
gauge-variant parts as 

(1) (1) (2) (2) (1) , 

(qa) =■■ Qa +£xqa, {qa) = -Qa +2£x (qa) + {£y - £x} qa, (2-36) 

(1) (1) (2) (2) ^^^ r 2 1 

(T^ab) =■ I^ab +£xT^ab, {'^ab) = '-nab +'2£x ij^ab) + {£y - £x]^ab- (2-37) 
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If we represent the multi-fluid system as an imperfect fluid system or we consider 
the micro-physical process, these energy flux and the anisotropic stress are related 
to the other fluid components. ^'^^^ In this case, the gauge transformations (j2-36p 
and (j2-37p should be derived from this relation. ^^-^ However, in this paper, we regard 
these energy flux and the anisotropic stress as fundamental variables for an imperfect 
fluid, phenomenologically. 

2.3.3. Scalar field 

Finally, we briefly summarize the perturbations of a scalar field (^ whose energy 
momentum tensor is given by 

fa ' = -g'^^a^Vc^ - ^5,' {g^^Vc^Vd^ + 2V{^)) , (2-38) 

where V{(p) is the potential of the scalar field (p. The fundamental variable of this 
system is the scalar field (p itself. In perturbation theory, we also expand this scalar 
field (p as 

^ = ip + \ipi + -X^ip2 + 0{\^), (2-39) 

where i-p is the background value of the scalar field (p. Further, as in the cases of 
the fluids, each order perturbations of the scalar field (p is decomposed into the 
gauge-invariant part and gauge-variant part as 

(^1 =: ipi + £x^, (2-40) 

^2=:V2 + 2£x0i + {£y-£\)v, (2-41) 

where 'pi and ip2 are the first-order and the second-order gauge-invariant perturba- 
tion of the scalar field. 

§3. Generic form of perturbations of energy momentum tensors and 

equations of motion 

Here, we consider the generic expression of the perturbations of the energy mo- 
mentum tensors and equations of motion for a perfect fluid ( ^3.ip . an imperfect fluid 
( N3.2p . and a single scalar field ( ^3.3p . We derive these perturbative expressions in 
terms of gauge-invariant variables defined in the last section. We also show that all 
perturbative variables are given in the same form as Eqs. (|2-20p and (j2-2ip . 

We note that we do not explicitly use any background values of the metric 
and matter fields in the derivations within this section. Further, we also note that 
we do not use any information of the Einstein equation nor the equations of state 
of the matter fields. Therefore, the ingredients of this section will be valid for 
any background spacetime and many perturbation theories of gravity with general 
covariance if the decomposition formula (|2-8p is correct. 

3.1. Perfect fluid 

The perturbative expressions of the energy momentum tensor for a perfect fluid 
are already discussed in KN2007.^^ Therefore, in this subsection, we just summarize 
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the definitions and results in KN2007 for the perturbations of the energy momentum 
tensor in ^3.1.11 In addition to the results in KN2007, we also show the perturbative 
expression of the equations of motion, i.e., the equation of continuity ( ^3.1.2p and the 
Euler equation ( NS.l.Sp . which are derived from the perturbations of the divergence 
of the energy momentum tensor. 

3.1.1. Perturbations of the energy momentum tensor 

The perturbative expansion of the energy momentum tensor (j2-22|) is given by 

(1) . (2) 

(p>f/ =: (p>r/ + A (p>r/ +-a2 (p>r/ +0{x^). (3-i) 

The background energy momentum tensor for a perfect fluid is given by 

ipyrj = {e + p)uau' + p6j. (3-2) 

(1) (2) 

The first- and the second-order perturbations ^^^T^ ^ and ^^^T^ ^ of the energy momen- 
tum tensor are also decomposed into the form as Eqs. (|2-20|) and (|2-2ip . respectively, 
i.e., 

(1) (1) 

(P>r/ =: (p)T„ ^ +£x^^^a^ (3-3) 

(2) (2) (1) 

{pYpJ =: (pyr^ ^ +2£x ^P^J + {£y-£x} ^^W (3-4) 

(1) (2) 

where gauge-invariant parts '^^ ^ and '^^ '' of the first- and the second-order per- 
turbations are given by 

(^) /(I) (i)\ (1) / (1) (1) \ 

(pyrJ:=U+v\ uav'+ V 5 J' + {e + p)\ua W -n'Wua+ Uau'y (3-5) 

J^^ h /(2) {2)\ /{!) (1)\ /(I) \ /(I) (1)\ (1) 



n" 



/(2) (1) \ 

+ ie + p)uJ U' -2n'"' Uc +2n'"ncdU^ - C'^'ud 

(1) /(I) , \ (2) , (2) 

+2 (e + p)Ua{ U^ -n^'uc ]+{e + p)Ua u^+ V K\ (3-6) 



where we defined 



(1) , (1) (2) , (2) 

W:= g"^ Ub, W:= g"^ Ub . (3-7) 



We also note that the fluid four-velocities Ua and Ua should satisfy the normal- 
ization conditions of the four- velocity 



rUaUa = g^'UaUb = -1. (3-8) 
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These normalization conditions yield 

(^) 1 

U" Ua = -HabU'^u', (3-9) 

(2) /(I) \ /(I) \ 1 

«" Ua = -gcb { W -W^'ud 1 ( W -n'^'Ua J + -CabU^uK (3-10) 

We have to emphasize that the perturbative expressions (j3-3p and (j3-4p are not 
definitions but the results which are derived from the definitions (|2-26p - (j2-28p . These 
are natural results from general formulae (l2-20p and ()2-2ip . However, these results 
imply that the framework developed in KN2003 and KN2005 does work in the case 
of the perturbations of the energy momentum tensor of a perfect fluid. 

3.1.2. Perturbations of the continuity equation 

Here, we consider the perturbations of the continuity equation for the perfect 
fluid which is derived from u^Vfe^^-'T'^^ = 0. This equation yields 

c(P):=rVae + (e- + p)^=0. (3-11) 

The energy density e and the pressure p are expanded as Eqs. (|2-23p . (j2-24p . respec- 
tively. Further, as shown in Appendix [XJ the four- velocity u"" and the expansion 6 
associated with the four- velocity Ua are expanded as Eq. ()A-46P and ()A-85P . Through 
these equations, we obtain the perturbative expansion of the continuity equation as 



C^P^ = Ci^) + a(i)cJ^^ + ^a2(2)c(^) + 0(A3) = 0, (3-12) 
where the continuity equation of each order is given as follows: 

Ci^) :=u"Vae + (e + p)0 = O, (3-13) 

(Dci^) := u'^Va ^e + (5) V,e + (^l^ + P ) 9+ ^9 {e + p) = 0, (3-14) 

(2)c(^) := u^^a ^f +2 (u^) V, ^6^ + (u^) V,6 

^ /(2) (2)\ (1) /(I) (1)\ (2) ^ ^ ^ ^ 

+e{ € + p]+2 e { e + p ]+ e {€ + p) = o. (3-15) 



The gauge- invariant variables for each order perturbations of e, p, u"", and 6 
are also given by Eqs. (f2^ - (l2^ . (fA48]l . f09l) . (fA^89l> - (lA^ . In terms of 
these gauge-invariant variables, the first- and the second-order perturbations (|3T4p 
and (I3T5P of the continuity equation are given in the gauge-invariant form. First, 
we derive the gauge-invariant expression of the first-order perturbation ()3-14p of 
the continuity equation (f3TT]) . Substituting Eqs. (p^ . (|X48l) . and (|X89]) into 
Eq. ()3-14p . we can decompose the first-order perturbation ^^'Cq into the gauge- 
invariant and gauge-variant parts as 

(i)ciP) = (i)cJP) + £^cJP), (3-16) 
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where 



(1)^^^ := u'^Va ? + (u'' -TC^^uA Vae + (^^ + v) 9 + {e + p) O • (3-17) 

We note that Eq. (j3-16p has the same form as Eq. (|2-20p . By virtue of the background 
continuity equation (|3-13p , the first-order perturbation (|3-14p of the continuity equa- 
tion is given in a gauge-invariant form: 

(1)^^) = 0. (3-18) 

Further, through the definitions (|2-26p - (j2-28p and the decomposition formulae ()A-48p . 
(|A-49p . ()A-89p . and (|A-9ip . the second-order perturbation (^^Cq^^ of the continuity 
equation defined by (|3-15p is decomposed into the form 

(2)^(P) ^ (2)^(P) ^ 2£x^''>c!p^ + {£y - £l) C^i\ (3-19) 

where 

/(2) {2)\ (2) /(l) \ (1) 

+9i£ +Vj + {e+p) + 2lW -W^b Va £ 

(1) /(I) (i)\ 
+20^+7'. (3-20) 



We also note that Eq. ()3-19p has the same form as Eq. (l2-2ip . Through the back- 
ground equation (j3-13p and the first-order perturbation ()3-14p of the continuity equa- 
tion, the second-order perturbation (J3-15P of the continuity equation is given in the 
gauge-invariant form: 

(2)cJ^) = 0. (3-21) 

Thus, we have obtained the gauge-invariant form of the first- and the second-order 
perturbations of the continuity equation for a perfect fluid through the lower order 
equations without any gauge fixing. 

3.1.3. Perturbations of the Euler equation 

Here, we consider the perturbations of the Euler equations. The component of 
Vfe'P^a'' = orthogonal to u"' gives the Euler equation 

C^P^ := (e + p)ab + r^Qbc^aP = 0. (3-22) 

Here, the three- metric (jab, the acceleration vector ab associated with the four- velocity 
Ua are defined by Eqs. ()A-3p . ()A-45p . and ()A-65P in Appendix [Al Substituting 
the perturbative expansions (l2^ . (12^21) . ([A^ . (lA4l) . and (|A47l) into the Euler 
equation (l3-22p . we obtain the expansion form of Eq. ()3-22p as 

C(P) =: cj^) + X(^)c^f^ + ^xmcj^P) + 0(a3) = 0. (3-23) 
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Then, we obtain the Euler equation of each order: 



CJ^P^ := (e +p)ab + g^'^qtc^aP = 0, (3-24) 



(i)C7(P) := {e + p) {a,) + {^^ + P j a, + g'^'q.^Va % 



(1) 
+5"" {qbc) VaP - h'^'^qbcVaP = 0, (3-25) 



(2) 
P 



^^^Ci'^ := (e + p) (al) +2 (^l' + ^P^) (al) + (^f + ?^) a, + ,-,,.V. 

(1) (1) (2) (1) 

+2^'*^ ((7be) Va P +5"'= (%c) VaP " 2h^''qt,,Va P 

-2/1'^^ (g^c) VaP + ( 2/i"'^/i/ - /'^^ gfc,V,p = 0. (3-26) 



As in the case of the continuity equation, the first- and the second-order pertur- 
bations (|3-25p and (|3-26p of the Euler equation should be given in the gauge-invariant 
form. First, we consider the gauge- invariant expression of the first-order perturba- 
tion (13^25)1 of the Euler equation. Substituting Eqs. ([23]), (|2^ . (lA^ . and (|A^54l) . 
we obtain the expression of (^-'C^^ as 

(1)^(P)^(1)C(P) + £^CW^ (3.27) 

where 

(^)C(^) ■.= {e+p) X + (^£+ P ) a, + 9^^q,,Va V 

+5"' Qbc VaP - g''''g''ndeqbcVaP, (3-28) 

The equation (j3-27p has the same form as Eq. (|2-20p . By virtue of the background 
Euler equation ()3-24p . the first-order perturbation ()3-25p of the Euler equation is 
given in a gauge-invariant form: 

(1)CJP) = 0. (3-29) 

Second, through the definitions (j2-26p . (|2-27p of gauge invariant variables, and the 
decomposition formulae ([MI), dHO]), (1X8]) . (lA^ . (lA^ . and (lA^ of gauge- 
invariant variables, the second-order perturbation ^"^'Cfj' defined by ()3-26p is decom- 
posed in the form 

(2)^b) ^ (2)c(p) ^ 2£x^^'^CP + (£y - £2^) Cj^\ (3-30) 

where 

I0-, M (2) /(2) (2)\ (2) (2) 

(2)C(P) - (g + p) ^, + / ^ + p j „^ + g--q,^Va V +5"^ Q6c V,p 

^ /(I) (1)\ (1) (1) (1) (1) 

-9'"^g'"Cdeqbc'^aP + 2( g +V] Ab +25"'= Qbc Va P -27^"^gfe,Va P 



-27^"'^ Qfee Vap + iW^g'^HdcqbcyaP- (3-31) 
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The equation (l3-30p has the same form as Eq. (12-21I) . Through the background 
equation (j3-24p and the first-order perturbation (j3-25p of the Euler equation, the 
second-order perturbation (j3-26p of the Euler equation is given in the gauge-invariant 
form: 

(2)ci^^ = 0. (3-32) 

Thus, we have obtained the gauge-invariant form of the first- and the second-order 
perturbations of the Euler equation for a perfect fluid without any gauge fixing. 

3.2. Imperfect fluid 

In this subsection, we consider the perturbative expressions of the energy mo- 
mentum tensor and the equations of motion for an imperfect fluid which includes 
important effects in the recent cosmology. We derive the gauge-invariant part of 
the perturbations of the energy-momentum tensor (j2-29p for an imperfect fluid in 
^3.2.11 Then, we derive the gauge-invariant continuity equation for an imperfect fluid 
^3.2.21 Further, we derive the generalized Navier-Stokes equation in ^3.2.31 which 
corresponds to the Euler equation for a perfect fluid. 

3.2.1. Energy momentum tensor 

To derive the perturbations of the energy momentum tensor for an imperfect 
fluid, it is convenient to consider the perturbation of the contravariant energy flux 

(f =:q- + \ (q-) +-X^ (q^) +0{\^). (3-33) 

(1) (2) 

The perturbations {q"") and {q"") are given by the same procedure as the derivations 
of the perturbations of u"" in Appendix IA.3I and these are decomposed into gauge- 
invariant and gauge-variant parts: 



(1) (1) 

(q-) = Q- -n''\b + £xq'', 




(3-34) 


(2) (2) . . (1) ^ (1) 

(q-) = Q- -q^r^ - 2n'' Q, +2?^'^^?^,,(z'^ + 2£x {q") +£Yq'' - 


f2 a 

- J^xq 


,(3-35) 


where we defined 






(1) , (1) (2) (2) 

Q'':=9'''Qb, 0.":= 9"'' Qc . 




(3-36) 


Further, it J' := g^'^Ttac is expanded as 






(1) -1 (2) 




(3-37) 



and the similar procedure to decompose the perturbations of u"" into the gauge- 
invariant and the gauge-variant parts yields 

(1) (1) 

(^a ) = -^a -"^ '^'^ac + £x'^a ; (3-38) 
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(2) (2) (1) 

(1) 

+2£x in,') +£Yn,' - £lnj', (3-39) 

where we defined 

(1) (1) (2) (2) 

nj':=g''nac, n,'':= g'" n,a . (3-40) 

The perturbative expansions of the traceless property [the last equation in Eqs. (|2-32p ] 
of the anisotropic stress it J' are given in gauge-invariant forms as 

(1) (2) (1) 

7r/ = 0, 77,"=?^"%ae, n,''=7TacC''' + 2n"'nca-2n"'n/TTac, (3-41) 

(1) (2) 

where we used vr^" ={nj')={nj')= 0. 

The orthogonal condition (j2-3ip of the energy flux ija to the four-velocity Ua is 
also expanded perturbatively through Eqs. (j2-34p and ()A-46p . The perturbations of 
the orthogonal condition (j2-3ip are given in the gauge-invariant form as follows: 

u^qa = 0, (3-42) 

(1) (1) 

n'^ Qa = -Qa W +qan'''nb, (3-43) 

U"" Qa = -Qa W' +2qan^' Uf, -2(7,7^"^7^,fen'' + qaC^'ub 

(1) (1) (1) , 

-2 QaW +2 Qa n^'W (3-44) 

Here, we have to emphasize that we did not fix any gauge choice. The perturbations 
of the orthogonal conditions ()2-3ip are also decomposed into the gauge-invariant 
and the gauge- variant parts as the formulae (j2-20p and ()2-2ip . Since the gauge- 
variant parts of the perturbations of Eq. (l2-3ip are given by the Lie derivative of its 
background value and the first-order perturbation of Eq. (|2-3ip . the perturbations of 
Eq. (|2-3ip are necessarily given in the gauge-invariant form through the lower order 
perturbations of Eq. ()2-3ip . 

Similarly, the orthogonal condition [the second equation in Eqs. (j2-32p ] of the 
anisotropic stress itab to the four-velocity n" are also perturbatively expanded through 
Eqs. (IA-46P and (I2-35P and these are given by 

U^TTab = 0, (3-45) 

(1) (1) 

U" Hab = -W nab + n^'UcTTab, (346) 

(2) (2) (1) 

u'' Hab = -W nab + 2W Uc nab - 2Wncdu'^nab + C^'u^nab 

(1) /(I) , \ 

-2 Hab W -H'^^'ub , (3-47) 
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where we have used the background orthogonal condition (J3-45P for the anisotropic 
stress Ua and its first-order perturbation. Eqs. (j3-46p and (j3-47p are gauge-invariant. 
This is due to the fact that the gauge-variant parts of the perturbations of u"'Ttab ai'e 
given by the Lie derivative of the lower order perturbations of vf^Ttab as Eqs. (|2-20p 
and (j2-2ip and these gauge- variant parts vanish as in the case of the perturbations 
of the orthogonal condition vf^qa = 0. 

These orthogonal conditions ()3-42p - (l3-47p for the perturbations of the energy 
flux and the anisotropic stress are necessary to specify the independent components of 
the gauge-invariant parts of the perturbations of the energy flux and the anisotropic 
stress. 

Now, we consider the first- and the second-order perturbations of the imper- 
fect part '*'T'^^ of the energy momentum tensor for an imperfect fluid. Through 
Eqs. (1225]), (|A46l) . (IMip . ([3^ . and dMZD, the perturbative expansion of the 
imperfect part (j2-33p of the energy momentum tensor is expanded as 



(1) . (2) 

(*>f/ = (^>r/ + A (^>r/ +-a2 (^>r/ +0{x^), 



(3-48) 



where 



(*>r/ = g.n'' + n.g'' + 7r,^ (3-49) 

(1) (1) (1) (1) (1) (1) 

(*>r/ = qa {u') + (qa) u' + Ua {q') + (ua) q'+ (vr,"), (3-50) 

(2) (2) (1) (1) (2) (2) (1) (1) (2) 
(^>r/ = qa {u") +2 {qa){u'') + (qa) u' + Ua (g") +2 {Ua){q') + (Ua) q' 

(2) 
+ i^a") ■ (3-51) 

Substituting (12^261) . (lA48l) . (12^361) . (I3^ . and (13^381) . into (13^501) . the first-order 
perturbation of the imperfect part of the energy momentum tensor for an imperfect 
fluid is decomposed as 

(1) (1) 

{iYpJ = (% * +£x^'^TJ, (3-52) 

where 

(1) (1) (1) (1) (1) (1) 

^'h-a '' = qaU'+Ua q'+ Qa u' + Ua O!' -2U(,(7,)7^*<^+ 77,'' -TTacH'' . (3-53) 

Further, through Eqs. (I2^ - (|2^ . (l3^3iD - (12^361 ). (l3^ . (R48]l . and (R49]l . the 
second-order perturbation (j3-5ip of the imperfect part of the energy momentum 
tensor for an imperfect fluid is decomposed as 

(2) (2) (1) 

{i)T^ b ^ (iyj-^ b ^2i;^ ii)T^ ' + {£y- £x) ^'K \ (3-54) 



where 



^(2) (2) (1) (2) 

(*>r„ '' = qaU' -2qan'" Uc +2qaH''^HdcU' - qaC'"Uc+ Qa u' 
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-UaQcC^^ - 2Uan''' Qc +2UarL^'"Hcdq'^ 

(2) (2) (2) (1) 

+Ua Q' + Ua g'+ nj> -2n'' n,a -TTacC^'' + 2n"'nJ'7Tac 

(1) (1) (1) (1)(1) (1) 

+2 QaU^ -2n^%, Qa +2 Q^Ua -2ri!"'q, Ua . (3-55) 

Note that Eqs. ([3^ and (f3^ have the same form as Eqs. (j2^ and (j2^ . 
respectively. 

Together with Eqs. (j3-3p and (j3-4|) for the perfect fluid, the total energy momen- 
tum tensor of each order is defined by 

(1) -I (2) 

f/ =: T/ + A T/ +-A2 T/ +0{\^), (3-56) 

(1) (1) (1) (2) (2) (2) 

J" b ._ (pYp b _|_ {i)ji 6 rp b.^_{pyp b _|_ (iYp b^ rp b.^_{p)j, b _|_ {i)j, b ^ (3.57) 

Further, the first- and the second-order perturbations of the total energy momentum 
tensor are also decomposed into the gauge-variant and gauge-invariant parts: 



(3-58) 
(3-59) 



Through Eqs. (j3-5p and (j3-53p . the gauge- invariant part 7^ of the first-order per- 
turbation of the total energy momentum tensor is given by 

(1) (1) (1) 

r/ := (p>r/ + (*)r/ (3-6o) 

(1) (i)\ (1) / (^^ (1) \ 

£ +V] UaU^+ V 5 J' + (e + p) n, U^ -n^<^UcUa+ Ua u^ 



(1) 


(1) 




'a 


= Tj+£xTa\ 




(2) 


(2) (1) 




Ta'- 


= Tj+2£xT,' + [£y- 


(1) 



(1) (1) (1) (1) 

+qa U' + Ua q'+ Qa v" + Ua Q' -2n(,(7,)?^''^ 

(1) 

+ nj -TTacH'''. (3-61) 

(2) 

On the other hand, through Eqs. (j3-6p . and (j3-55p . the gauge-invariant part 7^ ^ of 
the second-order perturbation of the total energy momentum tensor is given by 

(2) (2) (2) 

r/ := (p)r/ + (^>r/ (3-62) 

(2) (2)\ /(l) (1)\ /(^) \ /(l) {1)\ (1) 

g +V] UaU^ + 2{£+v\ua\U^ -n^W \+2{£ +v\UaU^ 
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/(2) (1) \ 

(1) (^^l . \ (2) (2) ^ 

+2 (e + p)Ua{ U^ -n^'uc \+{e + p)Ua n''+ V K- 

(2) (1) (2) 

+qa U^ -2qaV}" Uc +2qan^'^ndcu' - qaC^W^ Qa u^ 

(1) (2) (2) 

-Uaqcd"' - 2Uan^' Qc +2Uan^'ncdq'^ + UaQ^+ Ua q^ 

+ n^ -271^ Ilea -T^acC^ + 271^ H^ T^ac 
(1) (1) (1) {1)(1) (1) 

+2 QaU^ -2n^^uc Qa +2 Q^Ua -2n'"'qc Ua . (3-63) 

3.2.2. Perturbations of the continuity equation 

The equation of motion VbT^ ^ = is decomposed into the tangential and normal 
components to the four- velocity Ua- The tangential component of VfeT^'' = to Ua 
is given by 

u'^Vbe + (e + p)9 + (fat + V,,^'' + it'^^Bab = 0, (3-64) 

where we have used Eqs. (jA-ip . ()A-2p . and ()A-65p . Eq. (j3-64p is the continuity 
equation for an imperfect fluid. The normal components of V^T^^ = to the fluid 
four-velocity u"' is discussed in ^3.2.3i 

The first two terms in Eq. (j3-64p coincide with Cq defined by Eq. (|3-lip and 
the perturbative expression of these two terms are given in ^3.1.21 Then, we consider 
the perturbative expression of the remaining terms in Eq. (l3-64p : 

't(i) -ab ( pr - \_X7 p^ \ _i piiic-fea! 



C^O^ ■■= r {qa-ab + ^aqb) + t^Q^'^T^abBcd- (3-65) 

We expand this Cq as 

(1) . (2) 
Of = Cf + A Cf +'-X^ Cf +0(A3). (3-66) 

The first term g^'^qafib in Eq. (j3-65p is expanded through Eqs. (|2-34p . ()A-47p . and 
fAMl) . The last term g'^'t'^TtabBcd in Eq. (f3^65l) is expanded through Eqs. (IMHll . 

(1) 

()A-47p . and (|A-68p . Therefore, to derive the explicit expression of Cq , Cq , and 

(2) 

Cq , we have to derive the perturbative expression of the term Vaqb- As in the case 
of the perturbations of the tensor Aab which is derived in Appendix IA.21 we can 
derive the perturbative expansion of the tensor VaQb as follows: 

/ (1) \ 

Va% = Vaqb + A Va (qb) -H^b ' i^] Qc 
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+ \>? ( Va {%) -2H,,'[h] (g,) -H,,'[l] q, + 2h'''HaM [h] qA 
+0{\^). (3-67) 

Substituting the expansion formulae (l2^3iD . (12^351) . dMZD, (R^Tll . (U^SOll . and 

(1) (2) 

(fA^68l) . the perturbations cj\ ci*\ and C^*^ of C^^^ defined by Eq. p^ are given 
by 

Cf = q'ak + Vbq' + Tr'^'Bab, (3-68) 

Q*^ = q" (ab) + (%) a^ + V" (g,) -^"'i^ba ' [h] Qc - habQ^a' - habV'q' 

+7^'"' Bab + {TTab) B''' - 2hbd7r a' B'"' , (3-69) 

•^^A (2) (2) (2) 

Ci^^ = q' (db) + (%) a^ + V" (qb) -q'^a'U - VWb 

(2) (2) 

+ {7,ab) B-'' + vr^'^ i?,, -2lbd7T''''Bj - g'^'H^.'il] q, 

(1) (1) (1) (1) 

-2/1"'' ((7,) at - 2/i'^''V, (qb) +2h'''H^,' [h] q, + 2/i-/i/(7,afe 

+2/1-/1, V,g, + 25 V (vrab)5,, -4/i/7r^'' B^ -4/i/ (^,5) i?'"' 
+4/ife//i^'^7r^''S,rf + 2hachbd7r''''B'''. (3-70) 

Through Eqs. dl^D, (|2^ . (IMTD . (|X29l) . (lA^ . and (|A72]l . the first-order 
perturbation Cq of Cq is decomposed as 

(1) (1) 
C«=C«+£xC«, (3-71) 

where 

C« := gf- Ab + Qba' + V' Qb -g'^'H^b " i^] Qc - 'HabQ^a'^ - Hab'^'^q' 

+7r'^'' Bab + IJab B^' - 2nbd7Ta'B''''. (3-72) 

Similarly, through Eqs. i^, (l2a(l . (fM6]l . (IM7l> . (fA^ . (IX30]l . (jMl, dMS]), 

(2) 

(|A-72p . and (lA-74p . the second-order perturbation Cq is decomposed as 

(2) (2) (1) 

C7« = C« +2£x C7« + {£y - £\) C«, (3-73) 
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where 



(2) (2) (2) 

Cj*) := q' Ab +Qba' + V' Qb -q^'a'Cab - ^""q'Cab 

(2) (2) 

+ Hab B-' + vr^'^ B^ -2/:,rf7r«''S/ - g-'H^,^[C]q, 

(1) (1) (1) (1) 

+2g''''H'"^Habd [H] qc + 2g''' QaAb -2W'''' Q, a, - iW'qa A 

-2n'''Va Qb -2g'''H^b"[^] Qc +2n''''H^,'[n] qc + 2n'''nMaa'' 

(1) (1) (1) (1) 

+2n'''nabVcq' + ig'^'g'" HabB^d -mbeg'^'n'^' B^d -^g^'^U'^ U^b B^d 

+mbdn''f7r'''Bcf + 2n'"'n^''7TabBcd- (3-74) 

(1) (2) 

Here again, we have confirnied that the perturbations Cq and Cq defined by 
Eqs. p-69p and (l3-70p are decomposed into gauge-invariant and gauge-variant parts 
as Eqs. (j3-7ip and (j3-73p . and these have the same forms as Eqs. ()2-20p and ()2-2ip . 
respectively. 

Hence, through Eqs. (I3-I2p . ()3-16p . ()3-19p . ()3-66p . each order perturbation of 
continuity equation for an imperfect fluid is given by the gauge-invariant form as 

Ci^) + C« = 0, (3-75) 

(1) (1) (1) (1) 

(1) (1) 
= CJ^)+C«=0, (3-76) 

(2) (2) (2) (2) / (1) (1) \ 

Ci^) + C7« = Cj^) + C« +2£x C^^ + C« \+{£y- £l) {c!f^ + C^^ 

(2) (2) 

= CJP)+C«=0, (3-77) 

where we used the background continuity equation (I3-75P in the derivation of Eq. (I3-76P , 
and used the background continuity equation (|3-75p and its first-order perturbation 
(I3-76P in the derivation of Eq. (I3-77P . 

3.2.3. Perturbations of the generalized Navier-Stokes equation 

The normal component of V^T^ '' to Ua for an imperfect fluid (j2-29p is given by 



(e + p)ab + g/ [yap + u^Vbqa + Vfevf^"] + q" {qj + B^b) = 0. (3-78) 

This equation (I3-78P is an extension of the Euler equation ()3-22p for a perfect fluid 
to the imperfect fluid case. In this paper, we call Eq. (j3-78p as "^/le generalized 
Navier-Stokes equation"}^' 

As seen in the case of the Euler equation (j3-22p for a perfect fluid, the first 
two terms of the generalized Navier-Stokes equation (j3-78p is given by Cc which 
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is defined in Eq. ()3-22p . We call this Cc as "t/ie perfect parf of the generalized 
Navier-Stokes equation (|3-78p . Using this perfect part, the equation ()3-78p is given 
by 

CW + (7» = 0, (3-79) 

where Cc is defined by 

Cji'^ ■■= rUcd (u'\7bqa + g'^^bTtae) + g'^Qe {qbcO + Bbc) ■ (3-80) 



- (j\ 

We call this Cc as "i/ie imperfect parf of the generalized Navier-Stokes equation 
for an imperfect fluid. Since the perturbative expansion of the perfect part Cc is 
already given in ^3.1.31 we may concentrate on the perturbative expansion of the 
imperfect part (j3-80p . 

As in the cases of the imperfect part of the continuity equation for an imperfect 
fluid, the imperfect part of the Navier-Stokes equation is perturbatively expanded 
as 



(1) ^ (2) 

2' 



(7« =: C7« + A C^ +\\^ C7« +0(A3). (3-81) 



(1) (2) 



Here, C^ , C'j , and C^ are given by the perturbative expansion of Eq. (j3-80p 
through Eqs. (|2^ . (g^ZD, dA^l), (jA^ . fOTJl . (lA^ . ([i^85l) . and the pertur- 
bation of the tensor Vb^^ae- To derive the perturbative expression of VfoVfae, the 
perturbative expansion (j2-35p of the anisotropic stress and the perturbation ()A-25P 
of the connection C^^ are used. Further, through Eqs. ([MD, (pTOjl . (IM6]l . (IMTl) . 
dMD, dMD, (1X^291) . (TOOD . (IX72]1 . (1X74)1 . (IXMD . dMH), we can show that the 
first- and the second-order perturbations of the imperfect part C^ are also decom- 
posed into the gauge-invariant and the gauge-variant parts as 

(1) (1) 



c« 


= (f+£xCf^, 




(3-82) 


(2) 


(2) 


(1) 




ct^ 


= C« +£yC^:^ - 


-£^ci*) + 2£xC7W. 


(3-83) 



These have the same form as Eqs. ()2-20p and (j2-2ip . respectively. 

(1) (2) 

The derivations of the gauge-invariant part of Q and Q are straightforward. 
Therefore, we just summarize the results, here. The details of these derivations can 
be seen in the Appendix in Ref. 19). Through Eqs. (I3^2D - (|3^ . (13^301) . and (IMTD . 
each order perturbation of the generalized Navier-Stokes equation for an imperfect 
fiuid is summarized as follows. The background generalized Navier-Stokes equation 
is given by 



c(-) + c« 


= {e+p)ab+ g^'^qbc^aP 






+qbd [u^^cq" + ^c^"") + q' {qcbO + Bcb) 


(3-84) 




= 0. 


(3-85) 
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The first-order perturbation of tlie generalized Navier-Stokes equation is given by 

(1) (1) (1) (1) 

d^^ + C« = Ci'^ + C« +£x {CJ,'^ + C«) , (3-86) 

wliere 

^}\ ^]\ (1) /(i) (i)\ (1) 

Cr + Cb =ie + p)Ab + i£ +Vjab + g'^'qbc^a V 

(1) (1) 

+qM W Vcq'^ - qM'H''ueVcq'^+ Qm u'Vcq"^ - WUbdu'V cqa 

—'H"' qbd^'^T^ac+ Qbd '^c'^ '^ — qbd'Hce'^'^TT ^ 

(1) 

+qMVc n'" -qbdH^'f m TTf, - g''UbdH,''f[n\7:af 

(1) (1) (1) 

-n''^qeqbcO+ Q" qbcO + q^ Qbc + q^qbc 

(1) (1) 

+ Q" B,i, - n^'qeB^h + q" B^b ■ (3-87) 

Tlien, the first-order perturbation of the generahzed Navier-Stokes equation is given 
in the gauge-invariant form as 

(1) (1) 
C(^)+C«=0, (3-88) 

where we used the background equation (j3-85p of the generahzed Navier-Stokes equa- 
tion. The second-order perturbation of the generahzed Navier-Stokes equation is 
given by 

(2) (2) (2) (2) 



(1) (1) 

+2£x I C^^^ + C« \+{£y-£l] (C(^) + C« ) , (3-89) 



(2) /{2) (2)\ (2) (2) 



where 

(2) (2) 

Cb^ + C? = {e+p)Ab + [ £ + 'V\ab+ g'^'qbc^a V +5^' Qbc V aP 

. /(l) (1)\ (1) (1) (1) 

-g'"^g^^Cdeqbc^aP + 2[g+v]Ai, +2^"^ Q^, V^ V 



(1) (1) , 

-27e'qb,Va V -2n''' Qfe, V,p + 27^-^ V'Wrfe'ZbcVaP 

(2) 

+qbdu'Vc Q^ -qMUcH'''f[C]q} + 2qMUcn}gH^''[Hy 
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(2) (2) 

+qM U' ^cq"" - qMC'''ueV,q''+ Qm u'V cq'' - C''' quU^V ^qa 

(1) (1) (1) (1) 

-2qMUcH'''f[n] Qf +2qM U' V, Q" -2qM U, H'''^[n]qf 

(1) (1) 

-2qMn'%e^e Q'' +2qM'Hceu'H''''f[H]qf - 2qM'H'' Ue V^q'' 

+2qunfg7i''fuaVcq'' + 2 Qm u'V^ Q'' -2 Qm u,H'''f[n]qf 

(1) (1) ^ (1) ^ ., (1) 
+2 QmW Vcq^ - 2 Qm H'^UeVcq'^ - 27eUbdu''V c Qa 

(1) 
+27eUbdu'Hacf[n]qf - 27eUbd U' Vcqa + 2U''UbdH'''UeVcqa 

(1) 

-2n'"^ Qm u'Vcqa + 2nafnf\du'Vcq'' - C^UbdV'T^ac 

(2) (2) 

+ Qbd V.vr'^^ - qMCceV'T^'^^ + qbdVc H^' 

-qMH^'f^Tif, + 2g'^UbdnfgH,, <^[n]^f' - g-\dH,'^[£]naf 

(1) 

+2qMnfgH,'3[n]7r''f - 2H'''' Qm V^T^ac + 2n''UbdH'''V c^ae 

(1) 
-27eUbdy' {Hac) +2n'^qMH^''f[n]-Kf, + 27eUbdH,'^[H]'^af 

(1) (1) (1) (1) 

-2 Qm HceV'TT'^' + 2 Qm V^ 71'^" -2 Qm H^^'^injirf, 

(1) (1) 

-2g'''' Qm H,^f[H]^af - 2g6dHceV^ TI'^^ +2%d?^,p?^P'^Ve7r"'= 

+2qMn''m^pVy^ + 2g^^qMWH,J[n]7Tfe + 2qMn'''Hecf[ny^ 
-2qMH'''f[n] 77/e -25'^'^gMF/^[7^] 77,/ -£''qeqbc9 

(2) (2) (2) 

(1) , (1) 

-2W Qe 5b,0 + 2H'fnfeq'qMe - 2Wqe Qm 

(1) (1) (1) (1) (1) 

-2Wqeqbc +2 Q^Qm + 2 Q'' qM 9 

(1) (1) (2) 

+2q'QM0 + Q''B,h-C'''qeBcb 

(2) (1) (1) 

+q' B,b -2n'' Qe B,h - 2n''qe B^ 

(1) (1) 
+2 Q'^Bcb +2n'"^ndeq'B,b. (3-90) 

Through Eqs. (j3-85p and (j3-88p . the second-order perturbation of the generahzed 
Navier-Stokes equation is given in the gauge-invariant form as 

(2) (2) 

C(P) + cf )= 0. (3-91) 
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3.3. Scalar fluid 

Here, we consider the perturbations of the energy momentum tensor (I2-38P and 
the Klein-Gordon equation of the single scalar field ip. Since the perturbative expres- 
sion of the energy momentum tensor of the single scalar field is already discussed in 
KN2007, we just summarize the formulae for the energy momentum tensor in ^3.3. li 
After that, we consider the perturbations of the Klein-Gordon equation in N3.3.2[ 

3.3.1. Energy momentum tensor 

Through the perturbative expansions (|2-39|) and (|A-47p of the scalar field ip and 
the metric, the energy momentum tensor (|2-38p is also expanded as 

r/ = T/ + A«(T/)+iA2(2)(r/)+0(A3). (3-92) 

The background energy momentum tensor T^ is given by 

Tj:= Va^g'^Vcip - \5j' (VeV^VV + 2V {^)) . (3-93) 

Further, through the decompositions (p^ . (|2-10p . (|2-40p . and (|2-4ip of the pertur- 
bations of the metric and the scalar field, the perturbations of the energy momentum 
tensor ^^> (T^ '') and ^'^' (T^ ^) are also decomposed into the gauge- invariant and the 
gauge-variant parts as 



W(T/)=:(i)T/ + i:xr/, (3-94) 

(^HtJ) =:('>T^'' + 2£x^'HtJ) + {£y - £l)Tj, (3-95) 



where the gauge-invariant parts '^'7^ ^ and ^^^ '^ of the first and the second order 
are given by 

-\Sa' (v,(^VVi - V^ifW^^VdV + V,99iVV + 2'^i|^) ' (3-96) 



2\ 

1 



+2VaV3lVVl - 2V a^xU^"^ cV + V^^^aVV 



-VV^CdcVV + 2VcV7iVVi - 2VcV7iW'^'Vd^ 
+VcV72VV + 2(/^2^ + 2(/j?^-^j. (3-97) 

We also note that in these derivations, we did not used the homogeneous condition 
for the background scalar field (^ and the decompositions (j3-94p and (j3-95p coincide 
with the formulae (|2-20p and (j2-2ip . respectively. 
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3.3.2. Klein-Gordon equation 

Here, we consider the perturbation of the Klein-Gordon equation 

dV . 



Through the perturbative expansions (j2-39p and (j2-7p . the Klein-Gordon equation 
(I3-98P is expanded as 



C^K) ■■= ^"^a^ - ^^i^) = 0. (3-98) 



(1) 1 (2) 

CiK) =■■ C^K) + A C^K) +77 A2 C(x) +0(A3), (3-99) 



where 



CiK) ■■= VaVV - ||m = 0, (3-100) 



(2) 

CiK) ■■= V^VavSs - 2/i'^^V,Vb(^i - 2g^'H,,''[h] V ^^i + 2/i"^F,, ^ [/i] V c'f 



d'^v 2 d^y 



<^27T3tM-('^i)'773tM = 0. (3-102) 



Here, to derive the perturbative expansion of the kinetic terms of the Klein-Gordon 
equation (|3-98p . we use the connection ()A-25P between the covariant derivatives Vq 
and Va as in the case of the perturbative expansion of the tensor Aab in Appendix 

El 

(1) (2) 

The first- and the second-order perturbations Cik) ^-iid ^(-ftr) of the Klein-Gordon 
equation are also decomposed into the gauge-invariant and the gauge-variant parts. 
Through Eqs. dMD, (l2T0D . (f2^ . (l2^iTD . (fA^29l) . and (lA^ . the first- and the 
second-order perturbations ()3-10ip and ()3-102p of the Klein-Gordon equation are 
decomposed as 

(1) (1) 

CiK) ='■ CiK) +£xCiK), (3-103) 

(2) (2) (1) 

CiK) ='■ CiK) +'i£x CiK) + {£y - £x) CiK), (3-104) 

where 

CiK) ■■= v^v.y^i - H^-^[n\v,^ - n'^^'Va^bv - 'Pi-g^i'p), (3-105) 

+27i''''H,h"['H]Vcip - /:"^V„Vfc(^ + 27i''dn'^''VaVb^ - 2n''''VaVb^i 



^^^. ^ / ^2d^y 



9^27-^(95) - ((^i)'^T^((^). (3-106) 
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Here, we note that Eqs. (J3-103P and (j3-104p have the same form as Eqs. (l2-20p and 
(j2-2ip . respectively. 

By virtue of the background Klein-Gordon equation (j3-100p and the first-order 
perturbation (j3-10ip of the Klein-Gordon equation, the first- and the second-order 
perturbation of the Klein-Gordon equation are given in terms of gauge-invariant 
variables: 

(1) (2) 

C(K)=0, C^K)=0. (3-107) 

§4. Energy momentum tensors and equations of motion in cosmological 

situations 

Here, we derive the explicit components of the energy momentum tensors and 
the equations of motion for a perfect fluid, an imperfect fluid, and a scalar field. 
In the derivation of these equations, we do not use any information of the Einstein 
equations to guarantee the validity of the formulae to wide applications. 

4.1. Perfect fluid 

4.1.1. Energy momentum tensor 

As shown in ^3.1.H the gauge-invariant parts of the first- and the second- 
order perturbations of the energy momentum tensor for a perfect fluid are given 
by Eqs. ([33]) and (^. From Eqs. (1243]) . (R42]l . and (|XT3]l . the components of 
the first-order perturbation of the energy-momentum tensor of the single perfect 
fluid are given by 

(1) (1) (1) / n^ (1) {i)\ 

(p)T/ = -^, (pyr^'=-{e + p)iD'v+V-u'\, (4-1) 

(^) / (1) (i)\ (^) (1) 

(P)7:^ = (e + p)( A V + V.j, (pyr^^ = V5,'. (4-2) 

Further, from the components of the first- and the second-order perturbations of the 
fluid four- velocity which are given by Eqs. (|A-12p - ()A-14p . the components (j2-13p and 
(|2-15p of the first- and the second-order metric perturbations, the components of the 
second-order perturbation of the energy-momentum tensor of the single perfect fluid 
are given by 

, (^ (2) / (1) (1)\ / . (1) (1) (1)\ 

^Kj "^ = - £ -2{e + p) iOiV +Vi\ Id' V +V-iy'\, (4-3) 




[ (2) (2) (1) (1) {1)(1)\ 

+ {e+p)[Di V +Vi-2<p Di V -2 ^ H , (4-4) 
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(2) 



/(I) (1)^ 



(1) 



(1) (1)^ 



-2[£+V\\D'v+V- v' 



(2) 



(2) (2) 



(I) / ^,; (1) 



(1)^ 



+ (e + p) < -H' V -V' + u' -2 <P \D' V + V 



+2 i-2 ^^ 7*^ xA {d, %^ + v] - l;,M I , (4-5) 

(2) / n) (i)\ / (1) (1) (i)\ (2) 

(pV?: ^ = 2 (e + p) A V + V^\[ A V + V^' -v^ + V 5,'. (4-6) 



In Eqs. ()4-3p - (j4-6p . the vector- and the tensor-modes of the first-order perturbations 
are included into our consideration, which are ignored in KN2007. 

4.1.2. Perturbative continuity equation of a perfect fluid 

Now, we derive the explicit expression of the first- and the second-order per- 
turbations (j3-18p and (j3-2ip of the continuity equation for a perfect fluid in the 
cosmological situation. Before the derivation of the perturbations, we first show the 
background continuity equation (|3-13p in the cosmological situation, which is given 
by 



aCl^P^ = dr,e + zn{t + p) = 0. 



(4-7) 



Next, we consider the first-order perturbation (j3-17p of the continuity equation: 



i(i)cjp) = d^ f +3n re + vj + (e + p) M ^v -3d^ 1^^ j = o. 



(4-8) 



where we have used Eqs. ([2031) . (|XT2]l . (IXTSll . (1X^931) . (jMl, e = e(r?), and (FFp . 
On the other hand, the explicit expression of the second-order perturbation 
(j3-20p in the cosmological situation is given as follows. Through Eqs. (|2-13p . ()2-15p . 
(iX42]) - (IX44l) . (R^ - (R^ . e = e(r?), (g^l), and (jM]), the second-order perturba- 
tion (^Cq of the continuity equation is given by 



(9-, (v) r. (2) /(2) (2)\ / (2) ^ (2) 



{e + p)-E 



(p) 



0, (4-9) 



where 



"0 



(1) ^ (1) 

6 Ip Or^ Ip 



(1) (1) 

2 ^ + 4> 



A V 



(1) 



(1) 



+ 



D'^ ^'v' + V^ 



(1) 



D, 




dr. 



(1) w^ 

Dk V + Vk 






(e + p) 
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-2 D' V +V' -iy'\Di £ +2 ( 35„ 'F -Av\ ig +Vj. (4-10) 

4.1.3. Perturbation of the Euler equation 

Here, we consider the exphcit expressions of the perturbations ()3-29p and ()3-32p 
of the Euler equation for the perfect fluid. 

First, we consider the background Euler equation (j3-24p in the cosmological 
situation. In this case, the integral curve of the four- velocity n" = g°'^Ua, whose 
component is given by Eq. ()A-12p . is a geodesic. Then, we obtain 

Ofe = 0. (4-11) 

Through the background Euler equation (j3-24p . this yields 

%"V,p = 0. (4-12) 

This is supported by the homogeneous spatial distribution of the background pres- 
sure. Therefore, the background Euler equation (j3-24p is trivial due to the facts that 
the matter distribution is spatially homogeneous. 

Next, we consider the components of the gauge-invariant first-order perturbation 
(|3-29p of the Euler equation. Since the background fluid four-velocity Ua is tangent 

to geodesies, i.e., a^ = 0, we can easily see that the r/-component u^^^'C^ of (|3-28p is 
trivial due to the equations u^qbc = 0, Eqs. ()A-17p . ()A-59p . and (j4-12p . On the other 
hand, the spatial component (i-component) of Eq. (j3-28p gives the Euler equation 
for the perfect fluid: 

(I)C(^) = {e + p) I (a, + H) (d, ^l^ + v\ + A S| 

(1) / (1) (i)\ 

+Di V +d^P Di V + Vi = 0, (4-13) 

where we have used Eqs. (|2-13p . (|AT7p . and (|A-60p . This equation can be decom- 
posed into scalar- and vector-parts 

{e + p)\ {dr, + H) Di ^v + A ^S ]+D,V +dr,pDi ^v = 0, (4- 14) 



(1) (1) 

(e + p) {dr, + n) H +d^p Vi = 0. (4-15) 

Since the background value of the acceleration a^ vanishes, the r/-component of 
(2)cJP) in Eq. (IMTD is given by 

. ^ /(^) (1) (^A 

(2)C(P) = 2(i)ci^) u' -A V - VM = 0, (4-16) 

where we have used Eqs. ^EiW), (|A49]) . (jA^ . (jAW]) . g^c = 0, p = p(r/), and 
(|¥T3]1 . On the other hand, through Eq. (1203]) . ([A^20l) . (lA^ . ([A^ . and (|X62]) . 
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the spatial component of the second-order perturbation (l3-3ip of the Euler equation 
is given by 

(%?^ = {e + p) [{d, + H) (d, ^l^ + v\ + A S^l 

(2) ^ / (2) (2)\ i„) 

+Di V +d^P I A V + VJ - Sl^' 

= 0, (4-17) 

„(p) 
where ^^ is the collection of the quadratic terms of the linear-order perturbations 

defined by 

Sf> ■=-2<p Dilv -{e + p) ^ 



f^'^- (1) ('A f (1) / (1) W\ ] 

-2 (e + p) i.-' -A V - VM <^ A y -A A V + Vi H 

/(I) (i)\ f (1) . / (1) (^A / (1) (1)' 

-2 ^ + T' < I?^ <? +aJ D, V; + VJ + ?^ L>i ^ + Vi 



/ (1) (^A (1) 

-2 AV + Vi a^p. (4-18) 

As in the case of the first-order perturbations of the Euler equation, the equation 
(|4-17p is decomposed into the scalar- and the vector-parts as 

(e + p) {{d^ + n) A ^v +Di ^S] + A V +dr,pD, %^= DiA-^D^sf, (4-19) 



(e + p) {d^ + n) V. +d^p Vi= Sf ) - DiA-^D^sf\ (4-20) 

4.2. Imperfect fluid 

4.2.1. Energy momentum tensor 

Here, we consider the components of the each order perturbation of the energy 
momentum tensor for an imperfect fluid in the context of cosmological perturba- 
tions. In the context of the cosmological perturbations, the background spacetime 
is homogeneous and isotropic and the energy momentum tensor on the background 
spacetime is described by a perfect fluid. In other words, we can choose the four- 
velocity Ua of the dominant fluid so that the energy flux of the fluid vanishes, 

qa = 0. (4-21) 

Further, at the background level, we may neglect the anisotropic stress 

Tlab = 0. (4-22) 
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Then, the background energy momentum tensor is given by (l3-2p . Choosing the 
coordinate system as Eqs. ()A-12p . the components of the projection operator qj' : = 
9 Qac is given through Eq. (|A-16P and the background energy momentum tensor is 
given by Eq. (j3^ . 

Next, we consider the components of the gauge-invariant part of the first-order 
perturbation of the energy momentum tensor for imperfect fluid. In the situation 
where Eqs. (l4-2ip and ()4-22p are satisfied, the properties (I3-43P and (l3-44p give the 
ry-component of the first- and the second-order perturbations of the energy flux are 
given by 

(1) (2) (1) / ^s (1) (i)\ 

Qr,= 0, Qrt= -2 Qi D' V +V -y'Y (4-23) 

(1) a) (2) (2) 

where we defined Q,a='- a Qa, Qa='- a Qa, and used Eqs. (|2-13p and (|A-13p . In the 
same situation, Eq. (j3-46p and (|3-47p are given by 

(1) (1) (2) (2) (1) / (1) (1)\ 

77^^= =i7i^, nr,r,= 0, ^r,^= -2 flji [d^ V + V^ - v^ (4-24) 

(1) (1) (2) (2) 

through Eqs. (|2-13p and (|A-13p . where we defined tlah='- «^ ^ah and tlah='- ^ ^ah- 
Further, the traceless properties (j3-4ip of the gauge-invariant parts of the first- and 
the second-order perturbations of the anisotropic stress are given by 

(1) (2) (1) (1) 

7^'4i=0> i'tlii=1x^'flki. (4-25) 

The components of the gauge-invariant part (j3-6ip of the first-order perturbation 
of the total energy momentum tensor for an imperfect fiuid are summarized as 

(1) (1) 

(1) / ,l^ (1) {1)\ (1) (1) (1) (1) 

T^'=-{e + p)\D'v+V-i^'\-Q\ 7:^ = ^7.^+77.^ (4-26) 

(1) (1) (1). (1) 

where we have defined Q*:= 7*-' Qj and 77^-':= 7-^ 77^^. On the other hand, in 
the same situation, the components of the gauge-invariant part (j3-63p of the second- 
order perturbation of the total energy momentum tensor for an imperfect fiuid are 
summarized as 

(2) (2) / . (1) (1) 

T/ = - f -2(e + p) \d'v+V' 

(1) (1) (i)\ 

V +2 V* - i/M , (4-27) 
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(2) /(I) (1)\ / (1) (1)\ , f (2) (2) (1) / (1) (1)\ 1 



(?) {i)(P (1) / . (1) (i).\ 

+ Qi-2<P Qi +2 ilij L»J ?; + V^ , (4-28) 

(2). /(I) (lA /(I) . (1) (^)\ 

T^^ = 2{£+v\W-D'v-V\ 

i . (2) (2) (2) (1) / (1)\ 

+ (e + p) <^ -D* V -V' + i^' -2<p Id' V + v 

+2(-2^^^v5i)(z./i^^;-s)} 

- Q* -2 U +2 1^ Q* +2 x^'Q/ -2 77^* A V + Vfe - W , (4-29) 



(2) / (1) (1)\ / . (1) (1) (1)\ (2) . 

^V ( . (1) (1) m (1) / (1) (i)\ 

(2) (1) (1) (1) (1) 

+ n^^ +Aq^n^' -2x'"'nim, (4-30) 

(2) 

where we have used Eqs. (12031) . ([205]) . (|¥23]l . (R42ll ~ (|I^l4ll . and defined Q*:= 

(2) (2) (2) 

7*^' Qj- andi7/:=7^'=iT,fc. 

4.2.2. Perturbative continuity equation of an imperfect fluid 

Here, we derive the explicit expression of each order perturbation of the con- 
tinuity equation (j3-75p . (j3-76p . and (j3-77p for an imperfect fluid in terms of the 
components of the gauge-invariant variables. 

Since we consider the situation where Eqs. (j4-2ip and (j4-22p are satisfied, the 
background continuity equation coincides with Eq. (|4-7p for a perfect fluid. In the 
same situation, the gauge-invariant expression (j3-76p of the first-order perturbation 
of the continuity equation for an imperfect fluid is given by 

Cj^) + C« J = d, ? +3n ^S + v) +{e + p) U ^v -3d, %^ 

(1) 
+D^ Qj= 0, (4-31) 

where we have used Eqs. ()A-76P and (|4-23p . 
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Next, we consider the explicit expression of the second-order perturbation (I3-77P 
of the continuity equation for an imperfect fluid in the situation where Eqs. ()4-lip . 
(j4-2ip . and (|4-22p are satisfied. In this situation, the gauge-invariant imperfect part 
(j3-74p of the second-order perturbation of the continuity equation is given by 



(2) 



(2) (1) (1) 



(1) 
ik 



(1) 



(1) 



(1) 



(1)^ 



a C'' = D' Qi +4 ^ D' Qi -2 x*« A Qk +2d„ Qi D' v + V 



(1) 



+2Qi{D'{2 ^ 



(1) 



(1) (1) 



•^ 



(1) 



(1)^ 



2 {dr, + 2n) D' V + V 



(1) 



(1) 



+2 Aik D' D'' V + V'^ 






(!)■ 



+ f, X 



ik 



(4-32) 



where we have used Eqs. ([2031) . (I¥23D - (I¥25D . (IXMl) . dHO]), and (lA^76l) - (lA^79l) . 
Together with the perfect part ()4-9p of the continuity equation, the explicit form 
of the second-order perturbation (j3-77p of the continuity equation for an imperfect 
fluid is given by 



(2) (2) 





alc(^)+C, 



(2) 



(2) (2) 



(2) 



(2)' 



where we defined 
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/ .(1) « (1 
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(1) / (1) (1)^ 
-2A Qkh^ y^- X*' 
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« / . (1) (1)^ 
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(1) 



1 (i)\ / • (1) ^^' 



(1) (1)^ 
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(4-33) 



(4-34) 
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4.2.3. Perturbations of the generalized Navier-Stokes equation 

Here, we derive the exphcit expression of each order perturbation (j3-85p . (j3-88p . 
and (j3-9ip of the generahzed Navier-Stokes equation in terms of the components of 
the gauge- invariant variables. 

First, in the situation where Eqs. (l4-2ip and (j4-22p are satisfied, the back- 
ground generalized Navier-Stokes equation (I3-85P coincides with the background 
Euler equation (j3-24p for a perfect fluid. Therefore, we obtain Eq. (j4-lip . again 
due to Eqs. ()4-2ip and (j4-22p . Hence, the background generalized Navier-Stokes 
equation (j3-85p is trivial as mentioned above. 

Next, we consider the first-order perturbation (j3-88p of the generalized Navier- 
Stokes equation. As seen in igXIl through Eqs. (I¥TT]1 . (I¥23]l - (|¥25D . (|Xa6l) . (RTfl) . 
and ()A-76p . we can easily see that the jy-component of the first-order perturbation 
(j3-87p is trivial. On the other hand, in the same situation, the spatial component of 
Eq. (|3-87p is given by 



(1) 



(1) 



&) + df^ 



(e + p) [{dr, + n)\D,^v +Vi 



+ A 



(1) / (1) (i)\ (J) , (1) 

+A V +driP Di V +Vi]+{d^ + m) Q, +D^ n,k (4-35) 



0, 



(4-36) 



where we have used Eqs. ([2031) . g^SD, (li^2D . (lA42]l . (lAa6l) - (lA48]l . ([A^ . (|A76l) . 
and ()A-93p . The scalar-part of the generalized Navier-Stokes equation (j4-36p is given 
by 



(e + p) I {dr, + n) A ^v +Di i^H + Di V +dnpD, 



(1) 

V 



+DiA~^D^ 



(1) (1) ■ 



0. 
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Subtracting Eq. (j4-37p from Eq. (j4-36p . we obtain the vector-part of the first-order 
perturbation (j4-36p of the generalized Navier-Stokes equation: 



(1) 



(1) 



(1) 



(1) 



(e + p) (a, + n) Vi +dr,p v^ + {d^ + ah) Q^ +d' n, 



ik 



-DiA-^D^ 



(1) 



(1) 



(a^ + m) Qj +0" n. 



jk 



0. 



(4-38) 



Finally, we consider the explicit expression of the second-order perturbation 
of the generalized Navier-Stokes equation given by (j3-9ip . In the situation where 
Eqs. (|4Tip . ()4-2ip and (j4-22p are satisfied, the jy-component of this generalized 
Navier-Stokes equation (j3-90p is trivial due to the i-component of the first-order 
perturbation (l¥B6]l . Actually, through Eqs. (|i^ - (l¥25]) . (|A47ll ~ (|X2T]l . ([A^ - 
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(IA^62]> . ^AT7\\ - ^Al9\ . and (IXMl) . the r/-component of Eq. (f3^ is given by 



(2) (2) 

C(P) + Cf ) = 2 



(1) (1) 



^(1) . (1) (1)^ 



0. 



(4-39) 



On the other hand, the i-component of the second-order perturbation (l3-90p of the 
generahzed Navier-Stokes equation is given by 



(2) (2) r / (2)\ (2)1 (2) 
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where 
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(4-41) 



Equation ()4-40p is also decomposed into the scalar- and the vector-part in the similar 
form to Eqs. (j4-37p and ()4-38p with an additional source term — S"- * . 
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4.3. Scalar field 

In the inflationary scenario of the very early universe, scalar fields play important 
roles which drive the inflation itself and generate the seed of the density fluctuations 
through their quantum fluctuations. Keeping in our mind the applications of our 
framework of the second-order perturbation theory to this inflationary scenario, in 
this subsection, we summarize the explicit expressions of the perturbative the energy 
momentum tensor and the Klein-Gordon equation. 

4.3.1. Energy momentum tensor 

Here, we derive the explicit expression of the components of the energy momen- 
tum tensors (l3-93p . (I3-96J) . (I3-97J) for a single scalar field with the potential l^(</>) in 
gauge-invariant form. 

In the cosmological situation, we consider the homogeneous background field: 

^ = 99(7?). (4-42) 

Through the background metric (|2-12|) . the components of the background the energy 
momentum tensor ()3-93p are given by 

1 ^^ ^2 



T^'^ = -^{d,^Y-V{^l r/ = = r/, (4-43) 

Ti' = ^ {{dvVf - 2aV(^)} 7, '. (4-44) 



Through the components (j2T3p of the gauge-invariant part of the first-order metric 
perturbation, the components of the first-order perturbation of the energy-momentum 
tensor of a single scalar field are given by 

^% ' = -^ [d.,vdvVi- S {d:,^? + "'I^V'i} > (4-45) 



(4-46) 



^% ' = -^D.^idr,^, (4-47) 



^% ' = ^li' [d,^d,^^- $ {d,^f - a^l^^i} . (4-48) 

Finally, we consider the gauge-invariant part ()3-97p of the second-order pertur- 
bation of the energy momentum tensor for a scalar field. Through the components 
(|2-13p and (j2T5p of the gauge-invariant parts of the first- and the second-order metric 
perturbations and the homogeneous background condition ()4-42p . the components 
of the second-order perturbation of the energy-momentum tensor for a single scalar 
field are given by 



(2>r^ '? = -_] d^ipd^ip2 - (d^iff <P +a^V>2-g- - ^dr,^ ^ dr^ifi + A{dr,ipf{^] 



2 
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^'^ " = -^ Y^V i^i^^ - ^A"^! ^^A + 2 A 
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/(I) (2) 



+(a^99)2 4(^)2- u'ul - cp + (5^(^1)2 - A^^iI^Vi 
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(4-52) 



In the components (j4-49p - ()4-52p . the vector- and the tensor- modes of the first-order 
perturbations are included into our consideration, which are ignored in KN2007. 

4.3.2. Perturbative Klein-Gordon equations 

Here, we consider the explicit expression of the gauge-invariant part of each or- 
der perturbation (j3-100p and (J3-107P of the Klein-Gordon equation in the context of 
cosmological perturbations. Since the background field ip in cosmology is homoge- 
neous and the background metric is given by (|2-12p . the background Klein-Gordon 
equation (13-lOOp yields 



dV 

a^C^K) = dl^ + 2nd^^ + a^-r—i(p) = 0. 

(1) 



(4-53) 



On the other hand, the gauge-invariant part Ctx) of the first-order perturbation of 
(7(7^), which is defined by Eq. (j3-105p . is explicitly given by 



-2-M=0, (4-54) 






where we have used the background metric (|2-12p . the components (j2-13p of the 
gauge-invariant part of the linear-order metric perturbation, the homogeneous con- 
dition ()4-42p for the background field (/?, and the background Klein-Gordon equation 
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Finally, we consider the explicit expression of the gauge-invariant second-order 
perturbation of the Klein-Gordon equation in Eqs. (|3-107p . The gauge-invariant 

(2) 

expression of C(/^) for the Klein-Gordon equation is defined by (j3-106p . Through 
Eqs. (I2-13J) . (I2-15J) . (|4-42l) . and the background Klein-Gordon equation (|4-53p and its 

(2) 

first-order perturbation ()4-54p , the explicit expression of C(^k) is given by 

2 (2) / (2) (2)\ 

- a Cl^K) = drf^P2 + 2T-Cdr,ip2 - ^'^2 - \dr^ <P -^35^ ^ 1 dr^Lp 

(2) dV d'^V 

+2a2 <j _(^) + aV2^^(v') - ^(K) = 0, (4-55) 

where we defined 

. 2^^) d'^y ^ . 2, ^2Q'^y , X 

-4a ^ ^\-QZ^W)-a (931) -^(y) 

(1) (1) (1) 

+2 (a^ + 27i;) v' Di^i +4u' dr^Diifi - 2 x'^ DjD.ifi 

+2 <^ -2 <J dn ^ +6 <Ir a,, ip- -D* ( # + 1^ j i/j 

_2a2 y'V, ^(^). (4-56) 

§5. Summary and Discussions 

In summary, we have derived the explicit expressions of the second-order pertur- 
bations of the energy momentum tensor for a perfect fluid, an imperfect fluid, and 
a scalar field. Further, we also derived the explicit expression of the second-order 
perturbations the continuity equation and the Euler equation for a perfect fluid, 
the continuity equation and the generalized Navier-Stokes equation for an imperfect 
fluid, and the Klein-Gordon equation for a scalar field. As in KN2007,^^ we have 
again confirmed that the general formulation of the second-order perturbation the- 
ory developed in the papers KN2003^) and KN2005^''^ does work and it is applicable 
to the perturbations of the energy momentum tensors and the equations of motion 
for matter fields in the cosmological perturbation theory. In the derivations of these 
equations, we have again seen that the decomposition formulae (|2-2Up and (|2-2ip 
of the perturbations for the arbitrary fields play crucial roles in the gauge-invariant 
perturbation theory. Since the general relativistic higher-order perturbation theory 
requires the long algebraic calculations, in many cases, it is difficult to have con- 
fidence in the resulting long equations. In spite of this fact, we showed that all 
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perturbative variables are decomposed into the gauge-invariant and gauge-variant 
variables in the forms ()2-20p and (j2-2ip . This implies that the decomposition formu- 
lae ()2-20p and (|2-2ip are useful to check whether the resulting equations are correct 
or not. This is the main point of this paper. 

As mentioned in Introduction (5U, there are some attempts of the derivations 
of the perturbative expressions of the evolution equations of matter fields to the 
second-order.^^) For example, Noh and Hwang^^-* also summarized the formulae 
of the energy momentum tensor, equations of motion for an imperfect fluid in the 
cosmological situation up to second order without any gauge fixing. However, they 
implicitly imposed so-called "normal frame condition" in their formulae. In our 

(1) (1) (2) (2) 

notation, this normal frame condition corresponds to Dj v + Vi = Di v + Vi 
= 0. Due to this condition, their formulae are not equivalent to ours. Further, their 
formulae include gauge degree of freedom and the resulting formulae have some 
complicated forms. On the other hand, in this paper, all formulae are given in terms 
of gauge-invariant variables and there is no gauge ambiguity in these expressions. In 
this sense, the formulae in this paper are irreducible. We also have to emphasize that 
we did not fix any gauge when we derive any perturbative formulae in this paper. The 
gauge-invariant treatments of perturbations are equivalent to the complete gauge- 
fixing method. Therefore, we have realized the complete gauge-fixing without any 
gauge fixing. This is an advantage of the gauge-invariant perturbation theory in this 
paper. 

As emphasized in KN2007, the key point of our procedure is in the assumption 
which state that we already know the procedure to decompose the first-order metric 
perturbation into the gauge-invariant and variant parts. Mathematically speaking, 
this assumption is expressed by the statement that there is a vector field X"" which 
is constructed from some components of the linear-order metric perturbation so that 
its gauge transformation rule is given by the second equation in Eqs. Ii2-y\} and the 
linear-order metric perturbation is decomposed as Eq. li2-8\) . As shown in KN2007 
this assumption is correct at least in the cosmological perturbation case. However, 
even in the cosmological perturbation case, homogeneous modes of perturbations 
are excluded from our consideration because we assumed the existence of the Green 
functions A^^, (A + 2K)~^, and (Z\ -|- 3K)~^. These homogeneous modes of the 
cosmological perturbations will be some dynamical degrees of freedom in compact 
universes and we cannot say that these modes are unphysical. If we want to in- 
clude these homogeneous modes into our considerations, separate treatments will 
be required. Besides this detailed problem to extend the domain of functions for 
perturbations to includes homogeneous modes, the above assumption is correct on 
some background spacetime other than the cosmological background spacetime.^^-* 
Therefore, we propose a conjecture that the above assumption is correct on any 
background spacetime. To clarify whether this conjecture is true or not, non-local 
arguments on the background spacetime will be necessary, since gauge-invariant vari- 
ables are non-local one by their definitions. 

Even if the assumption is correct on any background spacetime, the other prob- 
lem is in the interpretations of the gauge-invariant variables. We have commented 



Perturbations of Matter Fields in the Second-order Gauge-invariant ... 39 



on the non-uniqueness in the definitions of the gauge-invariant variables in §2.2[ 
Tliis non-uniqueness in the definition of gauge-invariant variables also leads some 
ambiguities in the interpretations of gauge- invariant variables. Although the pre- 
cise interpretations of the gauge-invariant variables will be accomplished by the 
clarification of the relations between gauge-invariant variables and observables in 
experiments and observations, some geometrical interpretations are given through 
the explicit expressions of the gauge-invariant variables for the perturbations of ge- 
ometrical quantities like the acceleration a^, the expansion 9, the shear CTab, and the 
rotation ujah- For example, it is well-known that the expressions ()A-59p and (|A-60p 
of the first-order perturbation of the acceleration vector give the clear interpretation 

of the scalar-mode <P of the first-order metric perturbation. If u"' is a tangent to a 
geodesic even in the first-order perturbation, Eq. ()A-60P gives the equation of motion 
of a free-falling object in expanding universe and we can regard the scalar function 

(1) (2) 

^ as the Newton gravitational potential. Further, the similar interpretation for # 
is also possible through Eq. (JA-62J) though the quadratic terms of the linear-order 

. , (1) (2) 

perturbations in Eq. (|A-62p are difficult to interpret. Similarly, ^F and ^ contribute 

to the expansion of the first- and the second-order perturbations of the expansion as 

(1) (2) (1) 

in Eqs. (|A-94p and (|A-95p . respectively, and the vector- and tensor-modes Vi, fj , Xij, 

(2) 

and Xij contribute to the perturbations of the shear tensor through Eq. ()A-106P and 
()A-109p . As seen in these equations, the contributions of the second-order gauge- 
invariant variables c^ , if' , fj , and Xij to the gauge-invariant part of the second-order 
perturbations of the geometrical quantities 9, (Tab, i^ab are similar to the contribution 

of the first-order gauge-invariant variables ^ , <]/ , fj , and Xij of the first-order met- 
ric perturbation to the gauge-invariant part of the first-order perturbations of these 

(2) 

geometrical quantities. Therefore, we may regard the gauge-invariant variables <P , 

(2) (2) (2) 

^ , i^i, and Xij for the second-order metric perturbation are natural extensions of the 

. , , W (1) (1) (1) 

gauge-invariant variables <^ , ]p , i/i, and Xij for the first-order perturbations to the 

second-order one and the geometrical interpretations of these second-order variables 
will be similar to those of the first-order one. 

In all derivations of the perturbative expressions of the energy momentum ten- 
sors and the equations of motion in fj3l we did not use any information of the Einstein 
equations nor equations of state for the matter fields. This implies that the formu- 
lae in ^ are valid for the very wide class of perturbation theories of gravitational 
field. However, we have to emphasize that the equations of motion derived in §5] 
are not able to solve by themselves. These equations include metric perturbations. 
Therefore, to solve these equations, we have to use the Einstein equations. Further, 
we have treated the anisotropic stress and the energy flux in the case of an imper- 
fect fluid as independent variables of other perturbative variables. If we specify the 
micro-physics, these variables are related to the other perturbative variables. '^^'^'^•* 
Therefore, in general, it is meaningless to try to discuss the solutions to the equations 
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derived in f|4] by themselves. To discuss the solutions to them, we have to treat the 
Einstein equations and we have to specify the interactions between matter fields at 
micro- or macro-physical level. 

Finally, we can show that these equations of motion are not independent of the 
perturbative Einstein equations in the case where the spacetime is filled with a single 
matter field. ^'^^ This is a natural result because the Einstein equation includes the 
equations of motion for the matter field through the Bianchi identity at least in the 
case of the single matter content. Therefore, we may concentrate on the Einstein 
equations when we solve the system of the single matter field. On the other hand, 
when we consider the multi-fiuids or multi-fields system, we will have to specify 
the interactions between these matter fields and to treat the equations of motion 
derived in this paper, seriously. Therefore, in the realistic situations of cosmology, 
the formulae summarized in this paper will play crucial role. 
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Appendix A 

Perturbations of geometrical quantities 



The explicit expressions of the perturbative equations of motion for fiuids in ^ 
are given by the perturbative expansions of the fluid component e, p, Ua Qa, and itab 
as Eqs. ([2%D - (|?^ . (|2^3D - (|M5]) together with the perturbative expansion (|2^ 
of the metric perturbation. Further, to derive the perturbations of the equations of 
motion for matter flelds, it is convenient to consider the perturbation of the tensor 
VaUfe. From the components of the tensor V aUh, we define the acceleration, the 
expansion, the shear, and the rotation associated with the four- velocity u^. The 
acceleration Of, are defined by 

ab ■■= u^VaUb (A-1) 

and the expansion, the shear, and the rotation are defined as the trace, the traceless 
symmetric part, the antisymmetric part of the tensor 

Bab ■= qa%'^^cUd = ^aUb + Uattb, (A-2) 
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respectively. In this Appendix, we show the perturbations of the tensor field V aUb 
and the acceleration, expansion, shear, and rotation. All of the perturbations of 
these geometrical quantities are also decomposed into the gauge-invariant and the 
gauge- variant parts as Eqs. (j2-20p and (j2-2ip . 

A.l. Components of the gauge-invariant three-metric 

Before discussing the components of the perturbations of tensor Bab, we first 
derive the components of the perturbation of the three- metric qab, which is defined 

by 

qab ■=gab + UaUb (A- 3) 

and it is expanded as 

qbc = qbc + X [Qbc) +2^ (Qbc) +0{X^). (A-4) 

From the perturbative expansions (j2-7p . (j2-25p . each order perturbation of the three 
metric is summarized as 

Qab ■= 9ab + UaUb, (A-5) 

(1) (1) (1) 

{Qab) ■= Kb + Ua {Ub) + {Ua) Ub, (A-6) 

(2) (2) (1) (1) (2) 

{qab) ■= lab + Ua {ub) +2 {Ua){Ub) + {Ua) Ub- (A-7) 

Further, substituting the decompositions (j2-8p . ()2-10p of the first- and the second- 
order metric perturbations and the definitions of the gauge-invariant variables for 
the first- and the second-order perturbations of the fluid four- velocity in Eqs. (j2-26p 
and (j2-28p into Eqs. ()A-6P and ()A-7p . we can show that each order perturbation of 
qab is decomposed into the gauge-invariant and gauge-variant parts as 

(1) (1) 

{qbc) = Qbc +£xqbc, (A-8) 

(2) (2) (1) 

{qbc) = Qbc +2£x {qbc) +{£y-£x] qbc (A-9) 

(1) (2) 

Here, the gauge-invariant parts Qab and Qab are defined by 

(1) (1) (1) 

Qab ■='Hab + UaUb+ UaUb, (A-10) 

(2) (2) (2) {1)(1) 

Qab ■■= Cab + Ua Kb + Ua Ub + 2 UaUb . (A-11) 

To derive the components of the gauge-invariant part (jA-lOp and (jA-lip of the 
perturbative three-metric, the components of the gauge-invariant part of the fluid 
four-velocity are necessary. Components of the background and the first-order per- 
turbations of the four-velocity are summarized as 

Ua = -a{dr])a, (A-12) 
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(1) (1) / (1) m\ .(1) 

(2) (2) / (2) (2)\ . . (2) 

Ua = Ur, {dri)a + a A V + Vj {dx')a, D' Vi= 0, (A-14) 

(1) 

where the r/-component of the first-order perturbation Ua is determined by Eq. ()3-9p 

(2) ^ ^ 

and the 77-component of the gauge-invariant part Ua in Eq. (|A-14p is given by 

(2) f/(l)\' (2) / (1) (1) (1)W .(1) « W\1 ^, ^ 

Ur, = al{<p\ - <p -\DiV +Vi-Vi\ \d'v+V-v'\\ (A-15) 

through Eq. (13-lOp . where we have used the components (I2-13P and (J2-15P of the 
first- and the second-order metric perturbations. 

The components of the background three-metric qab defined by Eq. (lA-SP is given 
by 

qah = o?-hj{dx')a{dx^)b =■ a^lah- (A-16) 

Through Eqs. (j2T3p . ()AT2p . and ()AT3p . the components of the gauge-invariant part 
(lATOp of the first-order perturbation of the three- metric, which is defined by (lATOp . 
are given by 

(1) (1) (1) , /(i) (1) (i)\ 

Qr>r, = 0, Q^r^%^ ^ W " A ^ " Vi J , (A-17) 



Q^j = nij = a'i-2^ ^ij+ Xij j ■ (A- 18) 

Similarly, through Eqs. (12-15p . (|AT2p - ()A-15p . the components of the gauge-invariant 
part of the second-order perturbation of the three- metric, which is defined by Eq. (jATip . 
are summarized as 



(2) , / (1) (1) (1)\ / . (1) (1) (1)\ 

Q^r, = 2a^ Di V + Vi - W \D'v+V-v'y (A-19) 

(2) (2) ^ f(2) (2) (2) (1) / (1) {1)\ 1 

Qvn = Qvi= a^ IVi -Di V - Vi -2 <!> \Di V + V^\^ , (A-20) 

(2) , f (2) (2) / (1) (1)\ / (1) (1)\ 1 

Q,j =a^{-2^Ir 7,j.+ {1^ +2\Di V +Vi\\ Dj v + Vj ] } . (A-21) 



A. 2. Perturbation of the tensor Aab = "^aUb 

Here, we consider the perturbation of the tensor defined by 



Aab ■■= ^aUb. (A-22) 
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The covariant derivative Vq associated with the metric (jab on the physical spacetime 
is related to the covariant derivative Va associated with the background metric gab 
as 

Aab = VaUb - Cl^Uc, (A-23) 

where the connection C^^ is given by 

C\b = \f'' i^agdb + VbOda - Vdgab) (A-24) 

=: XH^,^[h] + Ia^ (i?,/ [/] - 2h^''Habd [h]) + OiX'), (A.25) 

where we defined 

Habc [t] := - {Vatcb + ^btca " ^ dab) , i^^ab ' [A ■= Q^^Habd [A (A-26) 

for any tensor tab of the second rank. Through Eqs. (|2-25|) and (|A-25p . the tensor 
Aab is expanded as 

( ^'^ \ 

Aab = ^aUb + A Va (ub) -H^a^lH Uc 

+^A2 I V„ (ul) -2H,, ^ [h] (ul) - [h,, ' [I] - 2h'''Habd [h]) u, I 
+0(A^) (A-27) 

(1) 1 (2) 

= :Aab + \Aab+-\^ Aab+0{\^). (A-28) 

Further, through Eqs. (|2-8p and (j2-10p . the first- and the second-order pertur- 
bations of the connection ()A-25P are given by 

H^, ' [h] = H^, ' [H] + VaVbX'^ + R^^, ^X^ (A-29) 

H,,^[l] - 2h'''Habd [h] = H^,'[C] - 2n'''Habd [H] 

+2£xH,b"[fA + VaV^y" + r^i^aeb" 

-£x (VaVfoX^ + R,,, ^X^) . (A-30) 

The components of the tensor H^f^ '^ \H\ are summarized in Appendix of KN2007 and 
the components of K^^^ '^ \C}\ can be derived from the components of H^^ '^ \H\ by the 
replacements 

W (2) (1) (2) (1) (2) (1) (2) 

^-><^, Ui^Ui, ^^W, Xij^Xij ■ (A-31) 

Then, through the last equation in Eqs. (j2-26p and (|A-29p . the first-order per- 
turbation of the tensor Aab in Eq. (|A-28p is decomposed as 

(1) (1) 

Aab = Aab +£xAab, (A-32) 
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where 

Aab :=Va Ub -UcH^,'[n]. (A-33) 

(1) ^ , (1) 

Aab in Eq. (|A-33p is the gauge-invariant part of the first-order perturbation Aab of the 

tensor Aab and we have verified that the first-order perturbation Aab is decomposed 
as Eq. (12^ . Similarly, through the last equation in Eqs. (12^261) . (I2^ . (1X291) . and 

(2) 

(|A-30p . the second-order perturbation 71^^ in Eq. (|A-28p of Aab is decomposed as 

(2) (2) (1) 

^afe = Aab +£YAab - £xAab + 2£x Aab, (A-34) 

where 

(2) (2) (1) 

Aab ■■= Va Ub -tic^6/[>C] + 2UcH'''H,ad [H] " 2i7,, ^ [?^] U^ . (A-35) 

(2) (2) 

Aab in Eq. (|A-34j) is the gauge-invariant part of the second-order perturbation Aab 

(2) 

of Aab and we have verified that the second-order perturbation Aab is decomposed 
as Eq. (IMTD . 

Now, we consider the components of gauge-invariant parts of the perturbations 
of the tensor Aab- First, we consider the components of the background value Aab- 
Through Eq. (1A-I2p . we obtain 

Aab = aHjij{dx'-)aidx^)b = aH'jab, (A-36) 

(1) 
where Ti := dr^a/a- The components of the gauge- invariant parts Aab of the first- 
order perturbation of the tensor Aab are summarized as 

(1) 
Ann = 0' (A-37) 

(1) / (1) (1) (1) (i)\ 

Ani = a[ dnDi v +dn Vi +Di <p +H W , (A-38) 

(1) / (1) (1) (i)\ 

Ar, = -an [Di V +Vi- W , (A-39) 

(1) f (1) (1) (1) / (1) (1) ^ (i)\ 

Aj = a I D,Dj V +Di Vj -D(i Uj^ - In <P +27i ^ +dn ^ j Jij 

+n Xij +ldn Xij] ■ (A-40) 

(2) 

Finally, the components of the gauge-invariant part Aab of the second-order pertur- 
bation of the tensor Aab are summarized as 

(2) r / (1) (i)\ (1) (1)1 /(I) . (1) w\ 

Ann = 2a{^n\DiV+V^]+n W + A ^ > ^' "^^ ^ " vM , (A-41) 
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f [^ A2) n (2) (2)\ „^ W(l) 

/(I) . (1) (i).\ / (1) (1) (1) 

+ v^ -D^ V - VM 2DiDj V +2Di Vj -2d^ ^ -fij 

,0 (1) on WM 
+5»7 Xij -2-D(i J^j) 1 ^ , 

^ (2) (2) (2)\ (1) / (1) (1)\ 

+dr, V^ + A ^ +7^ W - 2 <^ ( A ^ +W W J 

29,, 0/ 7ij - d^ Xij -2Dii Uj] j \ , 



45 



(A-42) 



Ar,i = ai dr,Di V 



I .(1) « «^ 



(A-43) 



(2) 

^^2 7 — ^ 



(2) (2) (2) / (2) (2) (2)\ 



Ir, (2) ^^ (2) 

+ 2^v Xij +^ Xii 

(1) r (1) / (1) 



+ 



W o W\ o (1) (1)1 

+3H ^ +2d^ ^ I 7ij - dr, Xij -2H Xij } 



on (') ^n (') 



(1) 



(1) 



(1)^ 



(1) 



(1)' 



-n-fij hli^' +n-fij D'^ V + vM pfc V + Vk 



(A-44) 



where we have used Eqs. (I2-I3p . (|2-15p . ()A-12p - ()A-14p . and the components of 

(1) 
^ab'^ [7^] and H^^^'^ [C] summarized in Appendix in KN2007. The components Arjrj, 

(1) (2) (2) 

Air), Arir), and Air) are also derived from the perturbations of the identity u Ag^^, = 0. 

A. 3. Perturbation of the acceleration ai, = u^'VcUt, 

Here, we consider the perturbations of the acceleration which is defined by 



ab := u^VcUb = u^Acb- 



(A-45) 



To obtain the perturbations of the acceleration, it is convenient to introduce 
perturbations of the contravariant four-velocity u"" which is expanded as 



(1) 1 , (2) 



(A-46) 



To obtain the relation between the perturbations n" and Ua of the fluid four-velocity, 
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we have to consider the perturbations of the inverse metric g which is given by 

-ab ^ gab _ ^f^ab ^ 1^2 (2h''^h^^ - l^^) + 0{X^). (A-47) 

Through Eqs. (j2-25p and ()A-47p . we obtain the relation between perturbations of 
u" and Ua of each order. Further, through Eqs. ([MD, OTO]) . (|2%]) . and (l2^ . 
the first- and the second-order perturbations of m° are also decomposed into the 
gauge-invariant and the gauge-variant parts as 

(1) (1) 

(n«) = U" -n^^uk + £xw", (A-48) 

(2) (2) (1) 

+2£x ( g"^ {ub) -h^^b j + £yu'' - £\v^, (A-49) 

where we used Eq. (l3-7p . These formulae are also given in KN2007 and we also 
note that Eqs. (|A-48P and ()A-49P have the same form as Eqs. (j2^ and ([2^ . 
respectively. 

Through the expansion formulae (|A-28p and (|A-46p . the acceleration Of, is also 
expanded as 

ab = afe + A (ab) +-A2 {ab) ^OiX^) (A-50) 

and we easily see that each order perturbation of the acceleration a^ is given by 

ab := u^'Aab = u^VaUb, (A-51) 

(1) (1) (1) 

(ab) := n'^ A^b + (n*^) At, (A-52) 

(2) (2) (1) (1) (2) 

(ab) := n" Aab +2 {u") Aab + (u"^) Aab- (A- 53) 

Substituting Eqs. (|A-32p and (|A-48p into Eq. (|A-52p . we easily see that the first- 
order perturbation of the acceleration is decomposed into gauge-invariant and gauge- 
variant parts as 

(1) (1) 
{ab) = Ab +£xab- (A-54) 

Further, substituting Eqs. (|X32]l . (lA^ . (lA48]l . and (|A49]l into Eq. (|A^53l) . we 
easily see that the second-order perturbation of the acceleration is decomposed into 
gauge-invariant and gauge-variant parts as 

(2) (2) (1) 

{ab) = Ab +2£x {ab) +{£y- £x} ab- (A-55) 

(1) (2) 

Here, we have defined the gauge-invariant parts Ab and Ab of the first- and the 
second-order perturbations of the accelerations by 

(1) (1) (1) 

A := ^i" Aab + W Aab - 'WWAab, (A- 56) 
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(2) (2) (2) (1) 

Ab:=u^ Aab + W Aab - C^'u.Aab - 2H''' U, Aab + 2n'"'ncdu''Aaj, 

(1) (1) (1) 

+2 AabW -2 Aab n'^'u,. (A-57) 

We note again that Eqs. (|A-54p and (IA-55P have the same forms as Eqs. (l2-20p and 
()2-2ip . respectively. 

Now, we consider the exphcit components of the gauge-invariant parts of the 
perturbations of the acceleration a^ associated with the fluid four- velocity Ua through 
Eqs. ()A-56P and (|A-57p . First, we consider the components of the background value 
of the acceleration associated with the four- velocity 

afe = u^Aab = aUu^-fab = 0, (A- 58) 

where we have used Eq. (|X36]l . Next, though Eqs. (f243]l . (fXa2]) . (lA43]l . (fA^Ml) - 

(1) 
()A-40p . and ()A-56p . the components of the gauge-invariant part Aa of the first-order 

perturbation of the acceleration are summarized as follows: 

(1) 

Ar^ = 0, (A-59) 

(1) (1) / (1) (i)\ 

A = A ^ + (^r, +n)[Di V +V^\. (A-60) 

Finally, we consider the components of the gauge-invariant part ()A-57p of the second- 
order perturbation of the acceleration. By the direct calculations through Eqs. p-13p . 
([245]) . (13^ . ([Aa2]) - ([Aa4l) . (|A^ - (|A^ . and (lA4T]) - ([A44l) . the components of 

(2) 

the gauge-invariant part Aa of the acceleration are given by 

(2) (1) /(I) . (1) m 

A^ = 2Aiiu' -D' V - VM , (A-61) 

(2) / (2) {2)\ (2) 

A, = {dr, +n)iDi V + Vi + A <? 

-2 ^ I 2Di ^ +{drj + n)iDi V +Vi\\ 

/(I) . (1) (1)\ / (1) (1) (1)\ 

+2 u^ -A V -Vn i Di V- -DjDi v -Dj Vi . (A-62) 

We can easily check the components ()A-59p ~ (|A-62p are consistent with the pertur- 
bation of the identity u'^ah = 0. 

A. 4. Perturbation of the tensor Bab '■= Aab + UaCib 

Here, we consider the perturbations of the tensor field Bab defined by Eq. ()A-2p . 
which is also written by 

Bab = Aab + Uattb (A-63) 
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through the tensor Aab defined by Eq. (IA-22J) . The tensor Bab is also decomposed 
into the trace part and the traceless part. Further, the traceless part of the tensor 
Bab is also classified into the symmetric part and the antisymmetric part: 

Bab = l^QabO + ^ab + U)ab, (A-64) 

9 := it'Bab = V,n^ (A-65) 

O-ab ■■= B(^ab) - -^qabO, (A-66) 

^ab-=Byab]- (A-67) 

6, aab, and Oab are called the expansion, the shear, and the rotation associated with 
the four- velocity Ua, respectively. 

The perturbative expansions of the tensor Bab is given by 

(1) 1 (2) 

Bab = Bab + X Bab +^X^ Bab +0{X^). (A-68) 

On the other hand, the perturbative expansion of the tensor Bab is also given by the 
direct expansion of Eq. ()A-63p . Comparing these perturbative expansions, we obtain 

Bab = Aab + Uattb, (A-69) 

(1) (1) (1) (1) 

Bab = Aab +Ua ttb + (Ua) ab, (A-70) 

(2) (2) (2) (2) (1) (1) 

Bab = Aab +Ua [ab] + (^a) 06 + 2 {Ua){ab) ■ (A-71) 

Further, substituting Eqs. ([A32]l . (jJ^ . and (f2^ into Eq. (jJ^ . we can decom- 

(1) 
pose the first-order perturbation Bab of Bab into gauge-invariant and gauge-variant 

parts as 

(1) (1) 

Bab = Bab +£xBab, (A-72) 

(1) 

where the gauge-invariant part Bab is defined by 

(1) (1) (1) (1) 

Bab ■= Aab +Ua Ab + Ua ab- (A- 73) 

Note that Eq. ()A-72p has the same form as the decomposition formula (j2-20p . Fur- 
ther, substituting Eqs. (f2^ . (f2^28l) . (fA^34l) . (jJ^ . fA55]) . the last equation in 

(2) 

Eqs. ()2-26p and ()2-28p into Eq. (IA-7ip . the second-order perturbation Bab of the 
tensor Bab is decomposed as 

(2) (2) (1) 

Bab = Bab +2£x Bab + {£y " £x) Bab, (A-74) 

(2) 

where the gauge-invariant part Bab is defined as 

(2) (2) (2) (2) {1)(1) 

Bab ■■= Aab +Ua Ab + Ua Ofc + 2 UaAb . (A- 75) 
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We also note that Eq. (IA-74P has the same form as Eq. (I2-21I) . 

(1) (2) 

Now, we consider the components of the gauge-invariant parts Bab and Bah of 
the first- and the second-order perturbations of the tensor field Bab- First, we note 
that the background value of the tensor field Bah is given by Eq. ()A-36p , (jA-SSh , and 



Bah = oHlah- 



(A-76) 



Through Eqs. (|A42l) . (|A43]) . (|A^58l) - (|A^ . and (lA73l) . the components of the 

(1) 
gauge-invariant part Bah of the first-order perturbation of the tensor Bah are sum- 
marized as 



(1) 
i3, 



vv 



(1) 
B. 



rji 



0, 

aU Vi -Di V -Vi 



(1) 
--B 



jjj) 



(A-77) 
(A-78) 



(1) 



(1) . (V 



(1 



aiDilDj V +Vj] - D(^i !>■) -[n ^ +2n ^ +dr, ^ ] jij 



(1) 



(1) 



(1) 



+^{dr, + 2n)iiY 



(A-79) 



Finally, through Eqs. (lAa2]) - (lA45]l . (|A4Tl) - (U44ll . (IX59]) - (IX62]) . ^KlEh . the com- 

(2) 

ponents of the gauge- invariant part Bah of the second-order perturbation of the tensor 
Bab are summarized as follows: 



(2) / (1) (1) (1)^ 

Br,r, = 2aH \ Di V +Vi - Vi 



' .(1) « «^ 

D' V +V' - u' 



(A-80) 



(2) 

Bni = a 



2 u^ -D^ 



(1) 






Di Id,^v+^ 



m (2) (2)\ (i)(i) 
+n[Vi-DiV-Vi\ -2n cpW 



(2) 

Bin = "■ 



2 u^ -D^ 



(1) 






Di I Di ''v + Vi 



(1) r, W 1^ (1) 

^(i ^j) -"V '^ 7ij + l^dr^ Xij 



(A-Sl) 



(1) a (^) 1-^ (1) 



(2) 



m (2) (2)\ (l)(l) 



/ (2) {2)\ (2) r / (2) (2)\ (2)1 

A I I?i ^ + VJ - Z)(, I'/) - <^ ?^ f 2 IP- + ^ j + a^ i^ }> -iij 



(A-82) 
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/W (I) (i)\ / (1) (1) (I) (I) 

+ j,k _j^k V _ yM 2Dfe IP- 7i,- - 4Z?(, 0^ 7j> + 2D^i Xj)k -Dk Xij 



H7^J z^fc^^' +H7^j D^ V + VM LOfc 'z;' + V,. 



(A-83) 



(1) (1) (1) 
We also note that the components Br/rj, -Bj^, Brji for the first-order perturbation 

(1) (2) (2) (2) (2) 

Bab and the components Brjr], Bir^, Brji for the second-order perturbation Bab are also 
derived from the perturbations of the properties 

U^Bab = U^Bba = 0. (A-84) 

A. 4.1. Expansion 9 

Now, we consider the perturbation of the expansion 9 defined by Eq. ()A-65p . 
Through Eqs. (|A-47p and ()A-68p . the expansion 9 is expanded as 

e = 9 + X^e +^X^''9 +0{X^), (A-85) 

where each order perturbations of 9 is given by 

= g^'Bab = g^'Aab = VaU\ (A-86) 

9 = g"' Bab -Babh''\ (A-87) 

^^9 = g"' Bab -2 Bab h^'' + Bab{2h'''K' - r'>). (A-88) 

,, . ,, , (1) 

Through Eqs. ()2-8p and (|A-72p . the first-order perturbation 9 is decomposed as 

(1) (1) 

9 =:e +£xO, (A-89) 

(1) (1) 
where the gauge-invariant part of is given by 

=: 5"' Bab -BabH'^'. (A-90) 

Similarly, through Eqs. (fMl) . (f2T0D . (IA-72p . and (|A-74j) . the second-order perturba- 

(2) 

tion is also decomposed as 

(2) (2) (1) 

9 := e +2£x 9 + {£y - £\] 0, (A-91) 
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(2) (2) 

where the gauge-invariant part of is given by 

(2) (2) (1) 

e ■■= 5"' Bab -BahC'''' - 2 Bab H'"' + IBabU'^'n,''. (A-92) 

Here again, we note that the decompositions (|A-89p and (IA-9ip of the perturbations 
of the expansion 9 into the gauge-invariant and the gauge-variant parts have the 
same forms as Eqs. (|2-20p and ()2-2ip . respectively. 

Now, we derive the exphcit expression of the perturbations of the expansion 9. 
First, we consider the background value of the expansion ()A-86p . From Eq. ()A-76p . 
the background value 9 of 9 is given by 

1 d„a 1 da . . ^ 

= S-T-C = 3^ = 3-— = 3H, (A-93 

a a^ adr 

where H is the Hubble parameter of the background universe and dr = adrj. Second, 

through Eqs. ([243]) . and ([A^ - ([AT9]) . the gauge-invariant part (|X90]) of the first- 

(1) 
order perturbation 9 is given by 

= -iA V -2,n <P -3dn ^ I . (A-94) 

Finally, through Eqs. (|2-13p . (|2-15p . and ()A-76p - (|A-83p . the gauge-invariant part 
(IA-92P of the second-order perturbation of the expansion is given by 



(2) 1 

= - 

a 



. (2) ^ (2) (2) ^ (1) /(I) (1)\ 

A V -3dr, ^ -Sn +65^ ^ i<P -2 ip \ 
(1) (1) /(i)\^ ^'^l (1) ^^^(1) 



( I. (1) ^^l\ { /(I) w\ ( (1) (1)' 

+ [d^ v +V^ \ l2Dk{^ - ^ \ + {2dr, + Sn) Dk V + Vk 

^^•l / (1) (1) (^) (1)^ 

+ x''' 2A i^k -2D,Dk V -2A Vk -dr, xL 



(A-95) 



A. 4. 2. Shear Uab 

The perturbations of the shear tensor aab defined by Eq. ()A-66P are given as 
follows. Through Eqs. (IA-4p . ()A-68p . and (|A-85p . the perturbative expansion of the 
shear tensor Uab is given by 

^ab =■■ (^ab + A {(Jab) +-X^ {(Jab) +0{X^), (A-96) 

where 

o-afe := B(^ab) - -^QabO, (A-97) 
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(1) (1) 1 (1) 1 (1) 

(aab) ■■= S(a6) -3 {qab) - -Qab 9 , (A-98) 

(2) (2) 1 (2) 2 (1) (1) 1 (2) 

(aab) := S(a6) -pab 6 "- {Qab) 6 "- (^afe) 0- (A-99) 

Through Eqs. (|A^ . (|A^ . and ([A^ . the first-order perturbation (jJ^ of 
the shear tensor is decomposed as 

(1) (1) 

{aab)=Sab +£x(Tab, (A-lOO) 

where 

(1) (1) 1 (1) 1 (1) 

^ab-=B^ab) -3 Qab - -Qab 9 ■ (A-101) 

Similarly, through Eqs. (lA^ . (KM . (lA^ . (TFTJll . (1X89]) . and (fA^ . the second- 
order perturbation ()A-99P of the shear tensor is decomposed as 

(2) (2) (1) 2 

0'ab = ^ab +'^£x (cTafe) +£YCrab - £x^ab, (A-102) 

where 

(2) (2) 1 (2) 1 (2) 2 (1) (1) 

^ab ■■= B^ab) -^qab -- Qab - - QabO • (A-103) 

Here again, Eqs. (jA-lOOp and ()AT02|1 show that the perturbations of the shear 
tensor are decomposed into the gauge-invariant and the gauge-variant parts in the 
same form as Eqs. (l2-20p and (l2-2ip . respectively. 

Now, we consider the explicit components of the shear tensor of each order 
perturbations. First, from Eqs. (|A-16p . (|A-76p . ()A-93p . and ()A-97p . the background 
shear tensor is given by 

aab = 0. (A- 104) 

Second, through Eqs. (lA46]l - (lA48]l . (IX77l) - fF79]l . and (lA^ . the components of 

(1) 
the gauge-invariant part Eab of the first-order perturbation of the shear tensor are 

given by 

(1) (1) (1) 

(1) if 1 \ (1) (1) (1) 1 (1) 1 

Eij = a I (DiDj - -7,,Z\j v +D^, V,) -D^, ^ +-d^ Xij \ ■ (A-106) 
Finally, through Eqs. (lA46ll - (lA^2T]l . (IA^ - (II^ . and (IX93ll - (IX95]l . the compo- 

(2) 

nents of the gauge-invariant part Uab of the second-order perturbation of the shear 
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tensor are summarized as 



(2) 



vv 



(2) 



0, 

(2) 



IT) 



s., 



Tjl' 



(2) 



m .(1) «.\ (1) 

(2) 



1 \ (2) (2) (2) 1 (2) 

DiDj - ^lij'^j V +Z)(i Vj) -D^i Uj) +-d^ Xij 
+ ^ ( 2D^i Vj) -dr, Xij I 



/ (1) (')' 
+2 D^ V + V'^ 



lk{iDj) - -lijDk 



(1) (ly 



/ (1) (1) 



1 ;^ fn «^S^ 






+2 I/'' I -2 ilkiiDj) - -lijDkj ^ +D(, Xj)k -^Dk Xij 



(1) 



+ 3^' 



Ik 



(1) / (1) (1)\ 1 (1) 

-jijDi Uk +jijDi iDk V + Vfc + -'yijdr, xik 



-likljl [A V --dr^ qr 



(A-107) 
(A-108) 



(A-109) 



A. 4. 3. Rotation uJab 

The perturbations of the rotation defined by Eq. ()A-67p are given by as fohows. 
The perturbative expansion of the rotation is directly derived from the perturbative 
expansion ()A-68P of the tensor field Bab and u)ab is expanded as 



and we have 



(1) 1 2 (2) Sn 

^ab =■■ ^ab + A UJab +-A UJab +0{X ), 



^ab ■= 0[ab], ^ab- = [B[ab])i ^ab- = {-t>lab]) ■ 



(A-110) 



(A-111) 



We can also define the gauge-invariant expression of the first- and the second-order 
perturbation of the rotation: 

(1) (1) (2) (2) (1) / 2 \ 

^ab = ^ab +£x^ab-, '^ab=^ab +2£x ^ab + {^Y - £ x) ^ab, (A-112) 
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where the gauge-invariant variables for the first- and the second-order perturbations 
of the rotation are given by 

(1) (1) (2) (2) 

^ab=l3[ab], ^ab=B[ab] ■ (A- 113) 

From Eq. (|A-76p , the background value of the rotation is given by 

^ab = an-f[ab] = 0. (A- 114) 

(1) 

Through Eqs. ()A-77p - ()A-79p . the components of the gauge-invariant part f2ab of the 
first-order perturbation of the rotation are given by 

(1) (1) (1) (1) (1) 

Further, through Eqs. (lA-80p - (lA-83p . the components of the gauge-invariant part of 
the second-order perturbation of the rotation are given by 

(2) 

n^y^ = 0, (A- 116) 

(2) (2) /(I) .(1) m / (i)\ 

Qir, = - ^r,i= 2a v3 -D3 V " V^ D[, V,] , (A- 117) 

(2) f (2) / , (1) (^A / (1) (1) 

+dnV[ilj]k]\- (A-118) 
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